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PASSIVE  SONAR  BEARING  ESTIMATION 
IN  THE  PRESENCE  OF  HIGHLY 
ANISOTROPIC  NOISE 
FIELDS 

By  L.  R.  Howell 

The  problem  of  sonar  bearing  estimation  in  the  presence  of  highly 
anisotropic  noise  fields  is  examined  for  both  conventional  and  optimal 
array  processors.  In  some  situations,  rhe  presence  of  one  or  more  (high 
level)  Interfering  noise  source  precludes  bearing  estimation  by  conven- 
tional, non-adaptive  array  processors.  However,  in  some  target  tracking 
problems  (e.g.  TMA)  the  performance  loss  due  to  noise  anisotropy  is  not 
total.  In  the  context  of  such  problems,  the  effects  of  multiple,  uni- 
directional interferers  on  the  conventional,  split-beam  tracker  and  the 
maximum  likelihood  estimator  are  investigated  for  linear  receiving  arrays. 

The  performance  of  the  maximum  likelihood  (ML)  estimator  is  charac- 
terized by  the  Cramer-Rao  variance  bound.  A new  result  is  obtained  for  this 
hound  which  is  valid  for  arbitrary  spatial  noise  distributions.  Expressions 
for  the  split-beam  tracker  (SBT)  output  mean  and  variance  are  derived  for 
the  noise  model  of  interest.  From  these  expressions,  approximate  theoret- 
ical measures  of  the  bearing  estimation  bias  and  standard  deviation  errors 
are  developed.  Further  statistical  measures  of  the  SBT  bias  and  standard 
deviation  errors  are  obtained  from  a digital  simulation  of  the  split-beam 
tracker.  The  effects  of  the  degree  of  noise  anisotropy  on  the  ML  and  SBT 
performance  measures  are  assessed  for  parametric  variations  in  the  aniso- 
tropic noise  component.  The  numerical  results  illustrate  the  complex 


behavior  of  random  and  bias  errors  with  variations  in  the  number,  separa- 
tion, and  position  of  interfering  noise  sources  when  such  interferers  are 
neither  coincident  nor  totally  resolvable.  As  might  be  expected,  anisotropic 
noise  components  at  remote  target-to-interf ercnce  separations  are  shown 
to  have  little  effect  on  either  estimator.  For  mean  separations  less  than 
about  two  resolution  bearowidths,  the  standard  deviation  of  error  for  the 
ML,  and  SBT  processors  is  primarily  dependent  on  the  spatial  distribution 
of  interference  power  in  the  vicinity  of  one  beamwidth  separatiop  from  the 
true  source  bearing.  The  relative  degradation  of  the  SBT  processor  over 
the  ML  estimator  is  similarly  dependent.  Multiple  interferers  bias  the 
SIJT  estimates  in  the  direction  of  the  mean  interference  bearing  with  the 
peak  bias  error  occurring  for  approximately  one  resolution  beamwidth 
separation  between  target  and  closest  interference.  Variations  in  signal- 
to-noise  ratio  and  array  size  are  also  examined.  Multiple  interfering 
noise  sources  do  not  affect  the  functional  SNR  dependence  of  the  selected 
performance  metrics.  A dramatic  decrease  in  SBT  bias  error  for  all  levels 
of  anisotropy  is  observed  as  array  size  increases.  The  standard  deviation 
of  error  also  decreases  with  array  size  with  the  rate  of  decrease  lower  for 
Larger  numbers  of  interferers. 
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INTRODUCTION 


] . 1 - Objective 

Due  to  its  crucial  role  in  "target"  localization,  bearing  estima- 
tion is  the  primary  post -detect ion,  sonar-array  processing  task.  The 
bearing  estimation  performance  of  optimum  and  sub-optimum  array  processors 
is  well  known  (see  references  [19],  [27],  and  [32])  for  the  case  of  a 
random  plane  wave  signal  corrupted  by  Gaussian  noise  which  is  isotropic 
and  uncorrelated.  In  many  applications,  however,  it  is  inevitable  that 
the  array  processor  will  encounter  anisotropic  noise  fields  containing  one 
or  more  directional  sources  of  interfering  noise.  It  is  the  purpose  of 
this  research  to  investigate  the  effects  of  a highly  anisotropic  noise 
field  on  bearing  estimation  by  optimum  and  sub-optimum  array  processors. 

A highly  anisotropic  noise  field  is  defined  as  containing  both  multiple, 
unidirectional  noise  sources  which  are  target-like  in  nature  and  an  addi- 


tive isotropic  noise  component  which  may  or  may  not  be  spatially  correlated. 


The  scope  of  this  study  is  limited  to  sonar  processor  structures 
which  are  of  a frozen  design,  i.e.  non-adaptive.  In  essence  we  are  assum- 
ing that  target  bearing  is  the  only  pertinent,  unknown  parameter.  This 
limitation  sidesteps  the  problem  of  defining  a suitable  methodology  for 
estimating  the  noise  field  parameters  (e.g.  the  Widrow  [35]  Griffiths  [11] 
or  Kelli  her  [10]  algorithms  for  optimum  processing  and  MacDonald's  [18]  sub- 
optimum  processors)  and  the  associated  problems  of  convergence  rates  and 
sensitivity  to  mismatch  of  assumed  and  actual  parameters  (see  Ccx  [7]  and 
Kooij  [17]). 
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1.2  - Background 

Early  studies  of  bearing  estimation  performance  [4],  [14],  [31] 


addressed  cases  of  signals  known  exactly  except  for  random  or  unknown  non- 
random  parameters.  Brennan  [4],  for  example,  assumed  a known  signal  with 
random  phase,  obtained  the  Cramer-Rao  bound  on  the  mean  square  error  of 
any  unbiased  estimator  and  compared  it  with  the  performance  of  a sub- 
optimal,  amplitude-comparison,  monopulse  processor.  Significantly,  he 
demonstrated  the  efficiency  of  this  sub-optimal  estimator  for  large  signal- 
to-nolse  ratios.  The  applicability  of  these  and  other  early  investigations 
Is  limited  to  environments  that  can  adequately  be  modeled  as  Isotropic, 
spatially  uncorrelated  Gaussian  noise.  In  a later  study,  Seldman  [27] 
examines  several  error  bounds  on  the  performance  of  known  signal-random 
phase  bearing  estimators.  In  addition  to  the  Cramer-Rao  bound,  Seldman 
obtains  the  Ziv-Zakai  bound  on  the  performance  of  any  arbitrary  estimator 
and  a modified  form  of  the  Ziv-Zakai  bound  which  offers  a larger  perform- 
ance bound  on  the  maximum- likelihood  class  of  estimators.  He  provides  an 
excellent  analysis  of  the  source  bearing  regions  over  which  these  bounds 
are  and  are  not  the  greatest  lower  bounds  and  those  regions  where  they 
are  Invalid.  Additionally,  the  low  SNR  threshold  below  which  the  Cramer- 
Rao  bound  is  inapplicable  is  presented  along  with  some  effects  that  array 
geometry  has  on  these  bounds.  However,  Seldman  also  restricts  his  analysis 
to  independent,  homogeneous  white  Gaussian  noise. 

MacDonald  and  Schulthelss  [19],  in  a comparative  analysis  of 
optimum  (maximum  likelihood)  and  sub-optimum  processors,  extend  the  bear- 
ing estimation  problem  to  broadband  signals  and  noise.  They  bound  the 
performance  of  the  maximum  likelihood  estimator  by  the  Cramer-Rao  bound 
and  derive  an  approximate  expression  for  the  rms  error  of  an  idealized 
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split-beam  tracker.  MacDonald  and  Schultheiss  also  Invoke  a spatially 
Incoherent  background  noise  model.  These  authors  emphasize  the  effects 
uf  array  size  and  geometry  on  the  sclit-beara  tracker  performance  relative 
to  the  C.ramer-Rao  bound.  Specifically,  it  is  shown  that  the  split-beam 
tracker  achieves  the  lower  bound  for  a two  hydrophone  array  while  for  a 
Large  number  of  uniformly  spaced  hydrophones  its  performance  is  degraded 
hy  approximately  4/3  . Furthermore,  the  performance  of  a multi-sensor 
array  whose  geometry  places  \ of  the  elements  at  each  end  of  the  array 
achieves  the  bound  exactly,  assuming  the  noise  remains  independent  at 
every  hydrophone. 

In  reality,  isotropic  noise  fields  exhibit  some  degree  of  spatial 
correlation,  not  only  for  special  array  geometries  such  as  MacDonald  and 
Schultheiss  [19]  suggest,  but  also  for  uniformly  spaced  arrays  whose  hydro- 
phone separation  is  small  relative  to  the  smallest  wavelength  of  Interest 
(e.g.  Jacobson  [15]).  Miller  [22]  investigates  the  effects  of  noise  corre- 
lation on  correlator  and  delay-9um  processors  for  the  case  of  known  signal. 
Miller  obtains  the  maximum- likelihood  bearing  estimates  from  the  outputs  of 
these  two  processors  and  concludes  that  the  delay-and-sum  processor  per- 
forms slightly  better  over  most  broadside  bearings.  Witt  [36]  also  exam- 
ines the  effects  of  spatial  correlation.  He  assumes  the  interelement  noise 
correlation  is  linear  and  extends  over  at  most  two  adjacent  hydrophones. 

The  performance  of  two  types  of  multiplicative  array  processors  la  analyzed 
from  several  theoretical  error  measures  and  from  the  results  of  twp  com- 
puter simulations.  Each  of  the  performance  measures  analyzed  Indicates  an 
Increase  in  standard  deviation  of  bearing  estimation  error  and  a afoift,  or 
bias,  in  the  direction  of  the  array  steering  angle.  Witt  concludes  that 
ihe  sum-and-correlate  estimator  is  generally  superior  to  the  monopulse 


i-st  inui:  <>r  over  the  range  of  noise  correlations  examinf5d.  Wirt  also 
i efereneos  some  additional  studies  which  compare  optimal  and/or  sub- 
optimal  process. >r  performances  in  isotropic  Dut  correlated  noise  fields. 

In  many  operational  situations,  it  is  impossible  to  describe  the 
noise  i toss-covariance  properties  at  the  hydrophone  array  as  simply  corre- 
lated and  isotropic.  For  example,  a noise  field  containing  many  unidirec- 
t lonal,  interfering  noise  sources  appropriately  describes  some  high  shipping 
density  areas  (e.g.  the  Mediterranean  Sea).  Also,  transitting  merchant  or 
military  vessels  present  a sonar  observer  with  multiple,  tightly  grouped 
interfering  noise  sources  of,  possibly,  similar  characteristics.  Similar 
highly  anisotropic  noise  models  are  appropriate  for  sonar  observation  of 


control ted-access  areas  such  as  harbors  and  channels.  One  can  also  en- 
vision several  geophysical  problems  wherein  multiple  directional  noise 
sources  are  likely  to  interfere  with  source  bearing  estimation. 

Such  highly  anisotropic  noise  fields  have  naturally  received  much 
attention  in  the  area  of  adaptive  array  processing.  In  general,  the  design 
of  adapt ive  array  processors  begins  with  limited  assumptions  of  the  signal 
and  noise  statistics  and  seeks  to  estimate  the  unknown  parameters  and  track 
any  subsequent  changes  in  them.  When  the  estimated  parameter  set  is 
.sufficiently  "close"  to  the  true  parameter  values,  one  is  assured  that  a 
chosen  estimation  criterion  (e.g.  maximum  likelihood  or  minimum-mean-square- 
error)  is  essentially  satisfied.  Questions  of  adaptive  array  performance 
commonly  address  the  rate  and  integrity  of  convergence  of  the  parameter 
estimates  and  the  sensitivity  to  mismatch  between  the  assumed  parameters 
and  their  actual  values. 

Static,  or  fixed  design,  array  processors  which  are  based  on  exact 
a priori  knowledge  of  signal  and  anisotropic  noise,  while  impossible  to 
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implement,  do  provide  Insight  into  the  asymptotic  performance  of  success- 
lolly  adaptive  implementations.  Hence,  there  is  a motivation  to  determine 
l lie  effects  of  highly  anisotropic  noise  on  static,  "optimum"  bearing  esti- 
mator structures.  On  the  other  hand,  sub-optimum  processors  are  often 
intentionally  designed  to  operate  in  isotropic  noise  environments  and  are 
commonly  of  a static,  or  fixed,  design.  Operational  situations  which 
violate  these  assumptions  of  isotropicity  motivate  analyses  of  anisotropic 
noise  effects  on  sub-optimal  processors.  In  spite  of  these  motivations, 
i lie re  are,  to  our  knowledge,  only  a few  analytical  studies  of  the  effects 
of  anisotropic  noise  on  the  bearing  estimation  performance  of  optimal  and 
sub-optimal  processors. 

Tn  a broad  study  of  the  maximum  likelihood  estimator  (ML)  and  two 
Interference-nulling  split-beam  trackers  (SBT)  MacDonald  [26]  conducts  a 
performance  analysis  of  noise  fields  containing  a single,  directional  inter- 
ference In  addition  to  isotropic  background  noise.  The  signal,  noise  and 
Interference  processes  are  defined  by  a set  of  non-random  parameters. 
MacDonald  presents  a lengthy,  theoretical  derivation  of  the  Cramer-Rao 
bound  on  ML  performance  and  approximate  expressions  of  bearing  estimation 
variance  for  one  pre-beamforming  and  one  post-beamforming,  interference- 
nulling, SBT  processor.  To  facilitate  analysis  of  these  somewhat  compli- 
cated expressions,  MacDonald  restricts  his  attention  to  cases  of  lpw  signal- 
to-noise  and  high  interf erence-to-noise  ratios.  He  shows  that  these  assump- 
t ions  lead  to  a strong  dependence  of  the  CR  bound  on  target-to-interf erence 
separation  and  number  of  hydrophones.  The  input-nulling  (pre-beamf orming) 
sub-op t Jinal  processor  can  effectively  eliminate  interference  effecps  except 
for  very  small  separations  (less  than  one  beamwidth)  where  it  loses  all 
capability  to  estimate  source  bearing.  The  output-nulling  (post-beamforming) 
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processor,  on  the  other  hand,  provides  satisfactory  performance  at  small 
separations  hut  performs  much  worse  than  the  input-nuller  at  separations 
greater  than  about  one-beamwidth , There  is  no  discussion  of  the  nature  of 
the  bias  error  from  these  two  sub-optimal  processors  which  would  quite 
possibly  be  significant  at  the  SNR  and  INR  levels  studied.  MacDonald 
recommends  that  the  problem  of  multiple,  directional  noise  sources  be 
addressed  in  future  investigations. 

Cox  [/]  treats  the  problem  of  angular  resolution  of  two  closely 
spaced  sources  by  conventional  and  optimal  array  processors.  Resolution 
probLems  seek  to  determine  the  angular  separation  at  which  an  observed 
phenomenon  is  known  to  result  from  two  sources  rather  than  one.  Cox  shows 
the  resolution  power  of  the  optimal  processor  to  be  about  three  times  better 
than  that  of  the  conventional  processor  for  the  classical  Rayleigh  resolu- 
tion Limit  and  high  SNR.  Cox  indicates  that  high  SNR  levels  are  required 
for  the  resolution  of  two  closely  spaced  interferers.  This  is  a reasonable 
result  since  one  would  expect  to  be  able  to  accurately  estimate  the  bearing 
of  a single,  low-level  source  long  before  being  able  to  resolve  two  closely 
spaced  sources  at  the  same  bearing. 

McCarty  121]  considers  a special  case  of  noise  anisotropy  which 
includes  multiple,  directional  noise  interferences.  In  the  context  of 
air-tralflc  control,  he  investigates  the  effects  of  multiple,  narrpw-band, 
zero-menu  Gaussian  interferences  on  maximum  likelihood  estimation.  This 
problem  allows  McCarty  to  describe  the  unknown  interference  parameters 
as  random  variables.  In  particular,  lie  assumes  that  the  interference 
bearings  are  uniformly  distributed  in  azimuth  over  the  interval 
around  broadside,  and  shows  the  resulting  spatial  noise  correlation  matrix 
to  be  a real,  Toeplitz  matrix  of  zero-order  Bessel  functions.  Hence,  he 


adJrusses  a special  ca;=e  of  noise  anisotropy  which  reduces  to  the  well 
known  "flat land"  or  cyli ndr leal ly  isotropic  noise  model  115),  [17).  McGarry 
derives  the  Cramer-Rao  bound  for  narrowband,  isotropic  but  spatially  corre- 
lated, white  Gaussian  noise  which  is  then  applicable  to  his  assumed  noise 
model.  Although  the  bound  requires  numerical  inversion  of  the  noise  corre- 
lation matrix,  it  is  greatly  simplified  for  this  special  case.  Tfre  maximum 
likelihood  estimation  performance  is  seen  to  rapidly  deteriorate  for  target 
hearings  away  from  broadside.  Even  at  broadside  bearings,  the  performance 
is  heavily  dependent  on  interf erence-to-signal  power.  It  is  shown  in 
chapter  2,  however,  that  errors  exist  in  McGarty's  expression  for  the  Cramer- 
Rao  bound . 

1.3-  Approach 

The  effects  of  highly  anisotropic  noise  fields  on  source  bearing 
estimation  by  optimal  (ML)  and  sub-optimal  (SBT)  array  processors  are 
assessed  from  the  behavior  of  several  performance  metrics  subject  to  para- 
metric variations  in  the  directional-interference  component  of  the  acoustic 
field  at  the  face  of  a hydrophone  array.  The  directional  interference  com- 
ponent of  the  total  noise  field  is  defined  by  three  parameters  — phe  number, 
mutual  separation  and  bearings  of  the  interferers.  The  detailed  behavior 
of  the  bearing  estimation  performance  metrics  is  examined  for  a fixed, 
baseline  set  of  signal  and  array  parameters  followed  by  a restricted  analy- 
sis for  extended  SNR  and  array  size  parameter  values. 

Chapter  2 discusses  optimal  or  maximum  likelihood  bearing  estima- 
tion and  its  performance  characterization  in  terms  of  lower  bounds  on  the 
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estimation  error.  The  Fisher-Dugue-Cramer-Ruo  (FDCR)!  bound  is  chosen  as 
an  appropriate  performance  measure  of  the  maximum  likelihood  estimator. 

A mathematical  model  of  the  signal-plus-noise-plus-interference  acoustic 
f ield-of -interest  at  the  hydrophone  array  is  then  presented.  Subsequently, 
a theoretical  expression  for  the  Cramer-Rao  bound  is  derived  as  a function 
of  the  inverse,  total-noise,  covariance  matrix. 

In  chapter  3,  the  split-beam  tracker  is  presented  as  a conventional 
sonar  processor  designed  for  maximum  performance  in  spatially  Incoherent 
noise  fields.  Given  the  signal-noise- interference  model  of  chapter  2, 
theoretical  expressions  for  the  mean  and  standard  deviation  of  the  processor 
output  are  derived  as  functions  of  the  array  steering  angle  and  the  method 
of  extracting  source  bearing  from  the  SBT  output  is  discussed.  This 
chapter  also  derives  an  approximate,  theoretical  expression  for  the  variance 
of  the  S1ST  bearing  estimate. 

As  is  often  the  case  with  analytical  studies,  it  is  difficult  or 
impossible  to  obtain,  under  the  necessary  controls,  pertinent  field  data 
for  comparison  with  the  approximate,  theoretical  results.  A successful, 
alternative  approach  to  field  experiments  is  the  exercise  of  a computer 
simulation  of  the  appropriate  array  processor  and  its  received  acoustic 
l lelds  (e.g.  Witt  [36]).  Chapter  4 describes  a digital  computer  simulation 
of  the  SBT  bearing  estimator  which  is  developed  to  provide  additional  meas- 
ures of  sub-optimal  processor  performance  in  highly  anisotropic  noise  fields. 
The  SBT  simulation  algorithm  is  discussed  along  with  the  methodology  employed 

'This  bound  is  usually  referred  to  as  the  Cramer-Rao  bound  al~ 
(hough,  as  pointed  out  by  Van  Trees  [33],  it  was  originally  formulated  by 
K.  A.  Fisher  and  proven  by  D.  Dugue.  With  due  respect  to  Mssrs  Flakier  and 
I'ugue,  we  too  shall  hereafter  refer  to  the  PDCR  bound  as  the  Cramer««Rao 
or  CR  bound. 
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tt>  generate  samples  of  the  input  signai-p  1 ••.s-noise-prus-interference 
piocesses.  An  algorithm  for  estimating  the  source  bearing  from  a zero- 
riossing  in  the  simulation  output  is  also  developed. 

l'he  effects  of  multiple,  directional  noise  sources  on  passive  bear- 
ing estimation  are  examined  in  chapters  5 and  b in  terms  of  the  parametric 
hi havior  of  the  Craaier-Rao  bound,  the  approximate,  theoretical  SBT  error 
mi  tries  and  the  statistical  measures  of  the  simulated  SBT  performances. 

The  total  directional  interference  is  defined  by  three  parameters  — 1)  the 
number  of  interferers,  2)  the  angular  separation  between  adjacent  Inter- 
t i rers  and  3)  the  mean  bearing  of  the  interference  cluster.  In  chapter  5, 
a fixed  baseline  set  of  signal  and  array  parameters  is  defined  and  the 
n suits  of  an  extensive  evaluation  of  the  performance  metrics  derived  in 
chapters  2,  3,  and  4 is  conducted.  In  chapter  6,  the  effects  of  extending 
tin*  baseline  values  of  SNR  and  array  size  are  examined  for  a more  restricted 
parametric  variation  of  the  total  interference  field. 

Chapter  7 summarizes  the  research,  presents  Its  conclusions  and 
oflers  recommendations  for  future  research  on  related  problems. 
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CHAPTER  II 


OPTIMUM  BEARING  ESTIMATION 


As  applied  to  parameter  estimation,  the  term  "optimum"  implies  that 
a properly  formulated  parameter  estimate  results  in  the  satisfaction  of 
some  pre-selected  performance  criteria.  For  example,  consider  Bayesian 
estimation  of  random  parameters.  Here,  one  selects  a cost  function  and 
defines  the  optimum  estimate  as  that  parameter  value  for  which  the  total 
risk  is  minimized.  In  source  bearing  estimation  problems,  we  are  interested 
in  unknown  but  non-random  parameters.  In  such  non-Bayesian  estimation  tasks, 
the  method  of  maximum  likelihood  (ML)  is  an  attractive  approach.  Very 
simply,  one  selects  as  the  ML  estimate  that  parameter  value  which  would 
most  likely  have  resulted  in  the  observed  data  set.  In  other  words,  the 
Ml.  estimate  is  the  parametric  location  of  the  peak  of  the  conditional  prob- 

-A  . 

ability  density  function  p(R  | b)  where  R is  the  observed  data  set  and  b is 
the  parameter-of-interest . Specifically,  the  ML  estimate  is  the  solution 
of  the  likelihood  equation 


tr  i^n  p(^  l b)i 


(2-1) 


which  occurs  at  the  largest  local  maxima  of  p(R  j b)  [33] . Equation  (2-1) 
is  similar  In  form  to  the  maximum  a posteriori  equation  in  Bayesian  esti- 
mation. Indeed,  the  maximum  a posteriori  probability  estimate  is  equal  to 
the  ML  estimate  when  the  amount  of  information  known  about  the  random  param- 
eter Is  negligible  (i.e.  infinite  variance;  completely  unknown  parameter). 
Consequently,  maximum  likelihood  estimation  has  an  intuitive  foundation  in 
Bayesian  estimation  theory. 
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The  appeal  of  maximum  likelihood  estimation  lies  largely  with 
several  asymptotic-  properties  that  it  possesses.  These  properties,  which 
apply  in  the  limit  as  the  number  of  independent  observations  increases, 
state  that  the  maximum  likelihood  estimate  is  1)  asymptotically  unbiased, 

2)  asymptotically  efficient  and  3)  asymptotically  Gaussian  distributed  [33], 
An  efficient  estimate  is  simply  one  whose  variance  achieves  the  Cramer-Rao 
lower  hound.  One  can  also  show  that  any  arbitrary,  unbiased  estimate  which 
is  efficient  is  identically  the  ML  estimate  [33].  Consequently,  it  seems 
reasonable  to  refer  to  maximum  likelihood  estimation  of  non-random  parameters 
as  "optimal"  estimation.  We  shall  subsequently  use  the  two  terms  inter- 
changeably. It  is  now  of  interest  to  address  the  performance  of  the  optimal 
estimator. 

2.1  - Tcrformance  of  the  Maximum  Likelihood  Estimator 

The  direct  approach  to  a performance  analysis  of  any  arbitrary  array 
processor  requires  that  one  first  obtain  theoretical  expressions  of  the  pro- 
cessor output  in  terms  of  the  assumed  input  signal  plus  noise  plus  inter- 
ference parameters.  One  can  then,  in  theory,  apply  a parameter  estimation 
algorithm  (e.g.  the  bearing  of  the  peak  processor  output)  and  obtain  ex- 
pressions for  the  estimation  bias  and  variance  errors.  This  direct  approach 
is  somewhat  naive  for  all  but  the  simplest  processors.  In  particular,  the 
solution  of  the  maximum  likelihood  equation  (2-1)  is  all  but  intractable. 

For  Ml.  estimation,  alternative  approaches  to  performance  analysis  have 
Included  use  of  simplifying  assumptions  and  theoretical  approximations  [13] 
[5],  numerical  solutions  of  equation  (2-1)  [9],  and  the  fortuitous,  a priori 
selection  and  analysis  of  a processor  structure  which  can  subsequently  be 
proven  efficient  [12].  However,  even  these  approaches  have  addressed  the 


more  tradable  isotropic  noise  probLera  and/or  small  (two  hydrophone) 
arrays . 

An  Indirect  approach  to  performance  analysis  of  bearing  estimators 
develops  mathematically  tractable  lower  bounds  on  the  estimation  error 
achievable  by  the  processor  of  interest.  This  approach  provides  a per- 
formance metric  which  can  be  evaluated  and,  of  course,  for  efficient 
estimators  it  reflects  the  actual  performance.  Two  primary  considerations 
govern  the  choice  of  an  error  bound  in  a processor’s  performance  in  any 
particular  problem.  The  first  consideration  involves  problem  dependent 
questions  such  as:  Is  the.  parameter  estimate  to-be-bounded  biased  or 
unbiased?  Will  the  estimation  error  be  confined  to  a localized  region  of 
the  parameter  space  or  must  global  errors  also  be  considered?  Secondly, 
one  obtains  the  best  indication  of  actual  processor  performance  by  studying 
the  behavior  of  the  greatest,  applicable  lower  bound  and  hence  seeks  to 
define  such  a measure. 

Unfortunately,  the  most  widely  applicable  error  bounds  do  not 
usually  prove  to  be  the  largest  lower  bounds  for  specific  problems,  In  a 
study  of  narrowband  signals  in  isotropic,  white  Gaussian  noise,  Seidroan  [27] 
evaluates  such  a bound  — the  Ziv-Zakai  bound.  He  defines  this  boynd  as  a 
lower  hound  on  the  average  mean  square  error  for  any  estimator  without 
restrictions  on  estimation  bias  or  ambiguities.  Seidman  also  presents  the 
will  known  Cramer-Rao  bound  and  a modification  of  the  Ziv-Zakai  bound  which 
Is  applicable  only  to  the  maximum  likelihood  class  of  estimators.  For  the 
hounds  Seidman  analyzes,  one  can  see  that  there  are  trade-offs  to  be  made 
prior  to  selecting  a particular  bound  as  the  performance  metric  of  optimal 
cst Imai ion.  For  example,  at  broadside  source  bearings  the  Ziv-Zakai  bound 
may  be  unnecessarily  conservative  while  for  bearings  approaching  endfire  the 


MI,  Ziv-Zakai  bound  becomes  Ineffective  and  the  CR  bound  asymptotically 
approaches  Infinity.  SLnce  the  CR  bound  applies  to  local  or  small-variara  c 
errors  it  is  also  shown  to  be  inapplicable  for  small  SNR.  Although  much 
Is  to  be  gained  computationally  from  analyzing  error  bounds  rather  than 
actual  processor  error,  it  is  apparent  from  Seidman's  work  that  ope  must 
carefully  select  the  most  appropriate  bounding  metric.  It  is  alsp  noted 
that  other  bounds  with  attractive  properties  (e.g.  the  Bhattacharyya  and 
Ikiranklu  bounds  [33]  and  the  J-divergence  [25],  an  information-theoretic 
distance  measure)  seem  almost  as  impractical  to  manipulate  as  true,  bearirg- 
estlmation  error  expressions. 

In  this  research,  the  bearing  estimation  performance  of  the  optimal 
processor  under  varying  degrees  of  anisotropic  noise  corruption  is  evaluated 
in  terms  of  the  parametric  behavior  of  the  CR  bound.  As  mentioned  above, 
the  advantages  of  employing  a bounding-type  metric,  in  general,  and  the  CR 
bound  in  particular  are  its  computational  simplicity  and  even  its  mathe- 
matical tractability  relative  to  both  other  bounding  metrics  and  the  actual 
Ml,  error  expression.  It  is  also  believed  that  the  CR  bound  is  especially 
reflective  of  the  ML  estimator  performance  for  the  problem  of  interest'  We 
are  primarily  interested  in  source  bearing  estimation  errors  arisipg  from 
i arget  tracking  situations.  Hence,  we  are  interested  in  true  source  bear- 
ings at  or  near  array  broadside  and  expect  to  avoid  the  asymptotic  growth 
of  the  CR  bound  near  endfire  bearings.  Furthermore,  since  we  have  chosen 
to  study  the  effects  of  anisotropic  noise  on  operational,  sub-optiipum 
processors  as  well  as  the  optimum  structure,  it  is  reasonable  to  restrict 
our  attention  to  low,  interf erence-to-noise  ratios  (INR=0.1)  and  moderate- 
lo-hlgh  s ignal-to-noise  ratios  (0.5^SNR<2.0) . We  have  some  assurance  in 
this  case  of  avoiding  problems  with  the  low-SNR  threshold  of  the  CR  bound 
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1 21 1 aiul,  according  to  Kl-Behrey  and  MacPhie  [9],  might  even  expect  the  Ml. 
estimate  to  bo  efficient  In  this  application.  We  shall  also  assume  in  this 
research  that  the  optimal  estimator  is  unbiased.  As  a final  motivation  for 
selecting  the  CR  bound  as  the  sole,  optimum-processor  performance  metric, 
we  reca  I L the  asymptotic  efficiency  and  unbiasedness  of  the  ML  estimate. 

We  can  thus  expect  the.  CR  bound  to  very  adequately  reflect  the  performance 
of  the  optimal  estimator. 

In  the  following  section,  we  develop  a mathematical  model  of  the 
acoustic  field  at  the  face  of  the  array.  In  section  2.3  we  derive  a new 
expression  for  the  CR  bound  on  the  bearing  estimation  variance  applicable 
to  arbitrary,  anisotropic  noise  fields. 

2.2  - Received  Signal,  Noise  and  Interference  Models 

In  target  tracking  applications  of  sonar  array  processors,  the 
received  acoustic  field  at  the  face  of  the  array  is  generally  assumed  to  be 
the  sum  of  a stationary,  Gaussian,  unidirectional  signal  process,  s(t,  0^.), 
and  a corrupting  (e.g.  self  noise,  ambient,  shipping),  Gaussian  noise  proc- 
ess, n (t),  which  is  possibly  spatially  anisotropic.  In  this  study  n (t)  is 
the  sum  of  M mutually  Independent,  unidirectional,  Gaussian  interferences, 

I (t,  0|n>,  and  an  independent,  isotropic  component,  n(t),  also  a Gaussian 
noise  process.  Hence,  the  total  field  at  the  face  of  the  array  is 

M 

x(t)  = s(t,  0T)  + n(t)  + ][  im(t,  em)  (2-2) 

m=l 

where  0 is  the  source  bearing  and  0 the  mth  interference  bearing.  Gen- 
T m — 

orally,  i distant,  interfering  noise  source  im(t,  9ffl)  will  be  non-stationary 
due  to  its  relative  motion  and/or  effects  of  the  propagating  medium.  Al- 
though In  certain  cases  it  is  feasible  to  mathematically  model  such  poise 
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processes  1 29 1 , It  is  beyond  the  scope  of  the  present  study.  Consequently, 
nil  ol  the  Caussian  processes  in  the  model  (2-2)  are  assumed  to  be  tem- 


poral I \ stationary. 

Several  simplifications  of  equation  (2-2)  arise  from  the  assumption 
(hat  all  directional  processes  originate  in  the  far  field  of  the  array. 

With  the  possible  exception  of  tow-ship  radiated  noise,  this  is  a quite 
reasonable  assumption  for  target  tracking  applications.  This  assumption 
allows  us  to  neglect  all  amplitude  gradients  across  the  array  respiting 
from  the  propagating  components  and  thus  reflect  the  directionality  of 
target  and  interferences  by  relative  propagation-time  delays  between  hydro- 
phones. Furthermore,  we  can  model  an  equi-phase  contour  of  the  propagating 
components  as  a plane  wavefront.  The  output  of  each  element  of  a hydro- 
phone array  Is  simply  a phase  delayed  version  of  equation  (2-2)  and,  for  the 
fill  sensor,  can  be  written  as 


x^(t)  = s(t  - A^)  + n^(t)  + 1 im(t  - s“)  (2-3) 

m=l 

where  is  the  source  delay  time  at  the  £th  phone  and  5^  is  the  mth  inter- 
i e renew  delay  time  at  the  £th_ phone,  both  relative  to  an  arbitrary  reference 
point  along  the  array.  Figure  2.1  Illustrates  the  acoustic  situation  mod- 
eled by  the  hydrophone  outputs  in  equation  (2-3)  for  the  case  of  a uni formly- 
spared,  linear  array.  Without  loss  of  generality,  we  can  define  hydrophone 
number  I as  the  phase  reference  and  obtain  the  following  phase  delays  at  the 
fill  hydrophone: 


CL^)d  . 

I c 

,m  a - l)d 
= • 

f c 


COS  0T 
COS  0 

m 


m * 1,  2,  . . . .M 


(2-4) 


where  d is  the  hydrophone  spacing,  c is  the  speed  of  sound  in  the  sea 

(r)000  ii/.;ec),  0 is  the  true  target  hearing,  M is  the  number  of  direc- 

i ional  interferences,  and  0 is  the  mth  interference  bearing. 

m — 

In  general,  the  isotropic  noise  component  at  the  £th^ hydrophone, 

( t ) , is  spatially  correlated.  Since  previous  investigations  [22],  [36] 
have  addressed  the  effects  of  isotropic  noise  correlation  and  since  our 
main  interest  Lies  in  the  effects  of  multiple  interferences,  n (t)  is 
assumed  to  be  spatially  uncorrelated  and  also  homogeneous  for  the  remainder 
of  this  research. 

The  signa l-plus-total-noise  field  output  from  each  hydrophone  is 
observed  by  either  optimal  or  sub-optimal  processor  only  during  the  finite 
time-interval  [0,  T] . Each  hydrophone  output  assumes  the  characterization 
of  equal  ion  (2-3)  over  the  T-second  observation  interval.2  If  one  allows 
i he  uniform,  linear  array  to  contain  L hydrophones,  then  the  total  array 
output  waveform  can  be  simply  expressed  in  vector  from  as: 

x(t)  = ^xj(t),  x2(t),  ....  x^(t)J  (0<t<T)  (2-5) 

Equations  (2-3),  (2-4),  and  (2-5)  along  with  the  aforementioned  simplifying 
assumptions  define  the  observed  signal-plus-total-noise  process  for  which 
the  optimal  and  sub-optimal  bearing  estimation  performance  metrics  will  be 
deve 1 oped . 


Although  we  shalJ  assume  the  hydrophone  outputs,  x^(t)  to  be  con- 
tinuous ior  all  t Lme  regardless  of  observation  interval  size  in  order  to 
tacilitate  Fourier  Series  representations,  Bangs  [2]  notes  that  precise, 

1 ternat  •-  representations  also  exist  when  x^(t)  is  assumed  to  vanish  outside 
t lie  [0,  T | observation  interval. 


I ii 


2.  3 - I iir  Cramcr-Rao  Bound*1 

ll  is  wi'Ll  known  [33],  [21] , [18j  that  the  Cramer-Rao  (CR)  lower 
hound,  ' , on  the  variance  of  any  unbiased  estimate  of  source  bearing  is 

"rv  n * <2’6) 


where  0 , is  the  unknown  source  bearing, 

\ 

x is  the  observation  vector, 

C(xjo^,)  is  the  log-likelihood  ratio  of  x given  6^  and, 

K [ • ] denotes  statistical  expectation  with  respect  to  the 
data  vector  x. 

Kvaluai ion  of  equation  (2-6)  first  requires  the  generation  of  ^(x|0^,)  which 


in  turn  is  defined  as 


<'(x|0T)  = log  A(x|0t)  = log 


jP(x j s + n, 

I p (x i N)  j 


(2-7) 


when-  | > ( :<  | S + N,  0^,)  is  the  probability  density  function  (pdf)  of  the  data 
given  1 1 1 a L a signal  at  bearing  0 occurs  in  addition  to  the  noise  yhile 
|)(x|N)  is  the  pdt  of  the  data  given  that  only  noise  occurs.  Since  the 
components  of  x (see  equation  (2-3))  are  sums  of  real  Gaussian,  random  noise 

—v  — V 

processes,  x wiJ I he  a real  vector  Gaussian  random  process  and  £(xj0^)  in 
(2-7)  loll ows  directly  upon  the  determination  of  the  mean  vector,  in,  and 
covariance  function,  K(t  - u),  of  the  data  vector. 

Generation  of  the  log-likelihood  ratio  is,  however,  facilitated  by 
a reformulation  of  the  data  vector  x(t)  into  a characterization  that  is 
orthogonal  In  the  temporal  dimension.  For  time-stationary  processes,  the 
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Fourier  .cries  of  x(t)  over  the  finite  interval  [0,  T]  is  such  a character- 
i/.atlon.  The  fill  component  of  x(t)  is  defined  as 


where 


A 1 . +jukt 

x(,(t)  = — - l X^(ja.k)  e 

k=-oo 

t -.it»kt 

X.,(jmk)  = / x^(t)e  dt 


(2-8a) 


k = ...  ,-1,0, +1, . ..  (2-8b> 


ire  the  Fourier  coefficients  of  x^(t)  and  where  ^ k is  the  Fourier 

conjugate  frequency  variable.  If  the  signal-plus-total-noise  process  is 
hand  1 lin i led  to  • m = 4-  — ^ N,  then  the  series  in  equation  (2-8a)  spans  only 
those  discreet  frequencies  (2N  + 1)  within  [-o^.+u^].  Furthermore,  since 
Xj,(i)  is  .i  real  function  of  time,  the  Fourier  coefficients  defined  by  (2-8b) 
are  conjugate  symmetric,  i.e. 


Xj,(ji»k)  - X^(-juk) 


k = 0,1, . . .N 


Consequenl ly,  (2-8a)  can  be  rewritten  as 

2 N ( +>kt 

x,(t)  = -y-  £ Re  jx£(>k)  e 

The  Fourier  characteri zation  of  the  observed  data  vector  can  now  be  ex- 
pressed as  the  following  N x L dimensional,  complex  vector  in  the  Fourier 
domain: 

I 

X^  = | (ja>o ),...,  X^(j^o)  » X\ (j  w j),..., X^(joi^  ),.. -Xj  (j  (2-9) 

where  Xp(j«)k)  is  given  by  (2-8b)  and  the  "prime"  superscript,  ’,  denotes  the 
vector  transpose. 

It  is  conceptually  straightforward  to  evaluate  the  CR  bound  for  the 
observation  vector  in  (2-9)  by  the  substitution  of  equations  (2-3)  and  (2-8b). 
The  end  result  Is  a theoretical  expression  which  reflects  the  dependence  of 
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lli.*  Cl!  bound  on  i lie  signal,  noise,  and  Interference  parameters.  Such  an 
ai'pro.uh  was  followed  by  MacDonald  1 18]  for  the  case  of  a single,  broadband, 
plane  wave  i liter  I crence.  The  algebraic  manipulations  required  to  effect 
tills  approach  for  a single  interference  are  quite  complex  and  lengthy.  Pre- 
liminary attempts  to  follow  this  approach  for  the  case  of  an  arbitrary  number 
ol  multiple  interferences  have  confirmed  MacDonald's  predictions  that  such 
a theoretical  derivation  is  essentially  intractable  (even  for  two  inter- 
lerences).  Consequently,  we  have  chosen  to  follow  McGarty’s  [21]  approach 
which  is  a mathematically  feasible  one  but  obscures  the  explicit  roles 
played  by  the  interference  parameters  since  it  requires  the  numerical  in- 
version of  the  noise-plus-interference  spatial  correlation  matrix  at  each 
and  every  frequency,  In  other  words,  the  isotropic  and  anisotropic 

noise  components  in  equation  (2-3)  are  lumped  together  into  a single  noise 

component  n „(t).  The  log-likelihood  ratio  and,  subsequently,  the  CR  bound 
a , 

are  derived  in  terms  of  the  inverse  of  the  covariance  matrix  of  n (t) . The 

a 

evaluation  of  the  bound,  discussed  in  chapters  5 and  6,  requires  the  numerical 
Inversion  of  this  total  noise  covariance  matrix.  Although  numerical  matrix 
inversion  is  a costly  operation.  It  will  be  seen  that  the  noise  covariance 
matrix  is  a complex,  Toeplitz  matrix  for  which  faster  inversion  algorithms 
(e.g.  1)0])  exist. 

The  derivation  of  the  log-likelihood  ratio  £(51^1 0,^)  for  Fovjrier 
characterizations  of  stationary  data  (e.g.  equation  (2-7))  is  straight- 
forward 1 21.  |,  [32],  Such  a derivation  using  our  notation  is  presented  in 
appendix  A for  convenience.  The  final  result  is  the  following  expression 
for  i lie  1 og-1  Ike  1 ihood  ratio: 


1 
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n c , •**  _i  j _i 

<n  "T>k-o\  ‘’l(u>k>  x lva(u,k)a  (v%  (2-10) 


-*  * 


f(X 


wlit?  I'  e 


°^(u,k)  = ” 


TS(aik) 


- 1 -v-  * 

1 + TS a (wk)a(odk) 

T is  the  total  observation  time, 

S(wk)  Is  the  signal  power  spectral  level  at  frequency  tok> 
X(wk>  is  a vector  of  the  L components  of  X,  (equation  (2-9)) 

corresponding  to  frequency  w , 

-1  k 

A“k>  - [F.{S(»k)B*(„k)}]  is  the  inverse  of  the  correlation 
matrix  (zero-mean  noise  processes)  of  the  Fourier 
characterization  (at  wk)  of  the  total-noise 
component  n *(t), 

tl)  V 

7?((nk)  is  the  phase  delay  vector  of  the  signal  process  with 

f tli  component  of  a(w  ) being, 

( ^ ) 

a !(  ^ = j " J wk  ~ cos  ®ijj»  £ = 1 » 2 » ••»»!*, 

(n ( • ) is  tlie  natural  logarithm,  and  * denotes  the  conjugate, 
transpose  vector  operation. 

Substitution  of  equation  (2-10)  into  equation  (2-6)  yields, 


CRA 


(2-11) 


where  C . (X(m  ) 1 0,^,)  is  the  quantity  in  braces  in  equation  (2-10).  Note  that 
0 enters  equation  (2-10)  only  through  the  signal  phase  delay  vector,  a(u^), 
and  always  appears  as  cos  (0  ) rather  than  9^.  If  we  let  p = cos  (9^,)  it 
can  he  shown  [5],  [20]  that  a relation  exists  between  the  variance  bound 
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.mi  (),„  ^ i von  by  (2-1])  and  a variance  bound  on  p(o^  j which  for  die  approx- 
1 CK,  p 

imately  broadside  bearings,  (0^,  c 90)  involved  In  target  tracking,  xs 


CR,  0 


CR,yi 


T sin"  9„ 


(2-12) 


Since  the  partial  derivatives  required  by  (2-11)  are  much  easier  to  obtain 
when  taken  with  respect  to  u rather  than  0T>  we  shall  develop  an  expression 

i or  o';.  and  Invoke  (2-12)  to  obtain  the  CR  bound  on  the  variance  of  source 

CR.ii 

bearings. 

We  can,  for  simplicity,  drop  the  explicit  reference  to  the  frequency 
dependence  of  the  component  terms  of  ^i(X(w^)  |0^,) , and  write  the  log- 
likelihood  ratio  for  the  ith  frequency  as 


\ • in(^s)+ 


l-A 

X 


(2-13) 


where  the  definitions  of  the  RHS  terms  of  (2-13)  follow  (2-10).  It  is 
interesting  to  note  that  the  likelihood  ratio  for  the  ith  frequency  given  by 
(2-13)  is  almost  identical  to  McGarty’s  [21]  result  for  narrowband  processes. 
II is  narrowband  result  differs  from  our  single-frequency  result  in  that  he 
assumes  the  number,  amplitude,  frequency,  and  bearing  of  each  inteyferer  to 
be  random  variables.  He  further  assumes  specific  pdf’s  for  these  parameters 
which  in  turn  result  in  a real,  Toeplitz  correlation  matrix  (K^  in  (2-13), 
however.  Is  complex)  comprised  of  zero-order,  modified  Bessel  functions. 

In  essence,  his  assumptions  define  a special  anisotropic  noise  model  equiv- 
alent to  the  well-known  "flatland"  or  cylindrically  isotropic  noise  model 


1151,  [17]. 
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I’lflorcM.i  iati"8  (2-13)  once  with  respect  to  p,  we  obtain 

i'"1  * S c" (s)+  T.'  **  *£'  * **  s^1  ■ 

+ rf~^  1 T?  "I 

T ^ X X *n  d + °Z  a X x K„  -r 


(2-14) 


XX%  ^ 


From  ih.  Jef ini t ion  of  a. 


(£)  - - 


,(£  - X)d  ( (P  - 'W a 

jo, exp  LL2. 


' ~ L)d  l 

\ 3W  c — 7 


_ , U - Dd 

- -jw  j—  <a)£ 


so  that 


<5a 


cl 


= (_jur7~)  ‘ dia8  (0,1,2, ...L-l) 


(2-15) 


= A a 


where  L li«-  matrix  A is 


A = (— jw  — — ) . 


diag  (0,1,2, .. .L-l) 


(2-16) 


A 1 so 


/ 


(A  a) 


?*)'■ 


— * A 

a A 


where  A is  the  conjugate  transpose  of  A defined  in  (2-16) 
properties  exist  for  the  matrix  A: 

A = -AC 
A = A' 

= -A* 


(2-17) 
The  following 


A 
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Differentiating  (2-14)  with  respect  to  p yields, 

■l  - - -2 


a _Li  3 2 , ( ’ t\ 

in7”  - 3^  f"  (tS) 


2 2 

0 °T  -k  * .)  ii  * 

+ 1T^  a **  xx 


■s1  • 


2 2 
daT  -i  * * _ 1 -v  _»-*  _ i _i  d°T  _j  * _]  ->  -x*  _i  _i 

+ 2 1‘ a A h x x h a + 2 -^a  h x x h * a 


-****_! 


-k  * _1  -I  * _i 


+oTa  AAK^  XXJ^  a + o‘  a XX  K,  AAa 

1 -1  -A*  _1  — k 

+2«2a  A XX  A a (2-18) 


Kvnluut ion  of  tlie  statistical  expectation  of  (2-18)  requires  knowledge  of 

l -x  _i  *. 

the  expectation  E 1 X X f . Substitution  of  (2-3)  and  (2-8b)  into  (2-9) 
yields 

_>  _k  M 

X = X (jw.)  = s(ju>  ) a + n (jw.)  + £ i (jw. ) b 
l l l in  l m 

m=l 

= s(joji)  a + ng  (ja>i)  (2-19) 

—X 

where  s<  j<n  ),  n(iii>j)  and  im  (ja>1)  are  Fourier  transforms  of  the  signal, 
noise, and  interference  processes  respectively, 

b|n  = exp  j-  j uj^  —■ — - • cos  0^  ^ is  the  phase  delay 

vector  of  the  mth  interference,  and 

.j 

n (jui)  is  the  Fourier  transform  vector  of  the  total 

fl 

corrupting  noise  process. 

Assuming  the  signal,  noise  and  interference  processes  are  band  limited  and 
assuming  T is  much  larger  than  the  correlation  time  of  these  processes,  the 
cxpectai ion  required  by  (2-18)  is  expressable  in  terms  of  the  power  spectral 
densities  of  the  component  processes  [3]  as: 
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r:  |X  X } = E j|i  a + n + [ [s  a + n + L i bm  J 

1 m=l  m=l  I 


= TS  a a “ + T N I + T I I b*  b 

— u m m • 


m=l 


m mm 


-i  * 

= TS  a a + Kjj 


(2-20) 


where  S,  N,  and  I are  Lhe  signal,  noise,  and  mth  interference  power  spectral 
dens (lies  at  the  implied  frequency,  and 

M 

K = T N 1 + T 7 I b b*  is  the  covariance  matrix  of  the 
■ N — *•  m m m 

total  noise  field  at  the  implied  frequency. 

The  matrix  is  Lhe  covariance  matrix  whose  inverse  must  be  computed  numer- 
ically in  evaluating  the  CR  bound. 

The  straightforward,  although  quite  tedious,  substitution  of  equation 
(2-20)  inlo  equation  (2-18),  yields 


(2-21) 


Details  of  the  above  derivation  are  carried  out  in  appendix  B.  Equation 
( 2-21)  is  an  expression  for  the  1th  additive  component  of  the  negative  in- 
verse of  the  CK  bound  and  is  valid  for  any  plane  wave,  unidirection3l , 

0. mss  Ian  signal  and  any  corrupting  Gaussian  noise  field  with  arbitrary 
spatial  power-distribution  defined  by  . Since  (2-21)  is  valid  for  any 
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arbitrary,  complex  covariance  matrix  iL  is,  to  our  knowledge,  an  original 

result.  McCarty's  comparable,  narrowband  result  [21,  eqn.  38]  which  is 

valid  only  for  real  inverse  covariance  matrices,  is  equivalent  to  the 

first  lour  terms  of  (2-21).  Not  only  is  equation  (2-21)  more  general 

than  his  result,  it  can  be  shown  that,  under  his  assumptions,  McCarty's 

-1  * 
result  is  In  error.  Since  his  K^,  is  a real,  symmetric  matrix  and  A = - A , 

McCarty's  result  requires  that  the  last  three  terras  of  our  equation  (2-21) 

vanish  as  a consequence  of  the  following  equality: 

A * + A = 0 (2-22) 


where  I)  denotes  the  null  matrix.  However,  observe  that 


* * • ”1  -1  * -1 
A ^ = (A  ) (J^  >'  « (K,,  A )’  = -(Kj,  A)' 


-1 

It  Is  easy  to  show  that  the  matrix  product  A is  symmetric  if  and  only 

- 1 


if  is  a diagonal  matrix.  Hence,  the  inequality  in  (2-22)  and  McCarty's 

expression  for  the  CR  bound  are  invalid  under  the  governing  assumptions. 
However,  if  the  sole  source  of  noise  is  the  uncorrelated  component,  the 
equality  (2-22)  would,  of  course,  hold.  In  cases  where  (2-22)  does  hold, 
not  only  do  the  last  three  terms  of  (2-21)  vanish,  but,  as  we  shall  see, 
the  first  two  terms  of  (2-21)  also  vanish. 

Kquatlon  (2-21)  still  requires  evaluation  of  the  derivatives  of 
Using  the  chain  rule  of  differentiation,  it  is  easy  to  obtain  the 
following  expression  for  the  first  derivative  of  a £ with  respect  po  p: 
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dy  _ / 

8 u 8'1 


(>~+  rTr*^) 

/ yv  y -»  * f * -1  -1  1 

--(o.p*a  [i  ^ 


(2-23) 


I ikewlse,  it  is  straightforward  to  obtain  the  second  derivative: 


fuT  <4  = -(a*)?j*A*  [a^1  + 2^a]  * + * *[aV  + J^a]  A a 

+ 2 (op3  ja  *[a*K^X  + K^1  a|  ~a  J (2-24) 


N.>w,  observe  tiiat  the  first  two  terms  in  (2-21)  can  be  written  as: 

,2  v x2_2 


a^ 


Jo2 

i a 0t  _ 

^4  alj7  (oT 


TS 

if 


/j,  °T  \ ^ °T  -v  * -1  t 

VnT S)+-^  » h * l 

l /^°tV  / -*  * -i  -i\  i 82°t 

if  M + r * *» 

T _ (84) 

8 7 \8y  / 


i | 82°- 

_ / 'PC  _ 


(<4>- 


(2-25) 


The  substitution  of  equations  (2-23),  (2-24),  and  (2-25)  into  equation  (2-21) 

yields. 


' jd7  Pi|  = "TS  j * *A*  [aV  + j^A  ] a + a * [a*K^  + K^a]  A a j 

+ 2TSai  !r  UV  + An  a]  <°t)2{"  %Xl+  S'a]  - f 

+ 2TSo2  JRe  A a)]  + |a  *A*K^1  t*\  2J 

- 2TS02  ja  [a  + K^a]  a j + of  1 1*A* [a*K^  1 + a 

-v  * r * _i  _i  i 

+ a [A  A*,  + A*,  A J A « 


(2-26) 


I <> ! low iii)',  final  expression  for  the  Cramer-Rao  bound  on  the  variance  of  any 
unbiased  estimator  of  source  bearing: 


win  re  each  variable  within  the  summation  is  dependent  on  the  implied  fre- 
<|u<  ncy  ut.,  and  where  G is  defined  by  equation  (2-27)  and  the  definitions 
<>l  l he  remaining  quantities  follow  equation  (2-10).  Recall  from  the  previous 
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<1  i sitiss  i .m  that  the  tractabil  ity  of  (2-30)  arises  from  our  willingness  to 
numerically  invert  the  noise  covariance  matrix,  K^,  at  each  of  the  N compo- 
nent f n quenc les.  The  explicit  definition  of  1°  terras  of  our  noise  model 

I o I lows  equation  (2-20)  above.  It  is  again  mentioned  that  "fast"  inversion 
algorithms  exist  which  exploit  the  Toeplitz  property  of  K^.  In  particular, 
when  numerical  inversion  of  is  required  in  this  investigation,  we  shall 
employ  the  algorithm  developed  by  Trench  [30]. 

MacDonald  [18]  was  able  to  derive  a theoretical  expression  for 
the  CR  hound  when  the  degree  of  anisotropicity  is  limited  to  a single  plane 
wave  ini  offerer.  MacDonald's  interest,  though,  is  in  the  effects  of  a 
single,  dominant  interference  on  a weak  source.  Since  our  result  in  equation 
(2-30)  i ; entirely  general,  it  is  of  interest  to  compare  its  behavior  with 
MacDonald's  theoretical  results.  Specifically,  he  presents  the  behavior 
of  a normalized  (.'R  bound  for  an  interference-to-noise  ratio  of  1.0  and 
I . ( S / N ) • • I,  where  I.  is  the  number  of  uniformly  spaced  hydrophones.  Figure 
2.2  compares  equation  (2-30)  with  MacDonald's  un-normalized  result  for 
variable  t a rgo t- t o- inter f erence  separations  and  for  the  parameter  values 
used  hi  llils  study.  Although,  for  the  parameters  of  interest,  there  is 
only  a minimal  hearing  estimation  capability,  figure  2.2  does  show  excellent 
agreement  between  the  two  expressions  for  the  CR  bound.  MacDonald  indicates 
that  In1  has  ignored  some  oscillatory  terms  in  his  graphical  presentation 
ol  (he  hound,  lie  assumes  these  terms  are  negligible  due  to  the  high  INR, 
low  SNR  condition.  Hence,  the  oscillatory  behavior  of  our  result  (equation 
(2-30))  in  figure  2.2  is  to  be  expected. 


CHAPTER  III 


SUB-OPTIMUM  ESTIMATION:  THE  SPLIT-BEAM  TRACKER 

In  general,  one  is  attracted  to  sub-optimum  implementations  of  para- 
meter estimators  by  their  potential  economic  benefits.  Of  course  such 
Implementations  are  usually  realized  at  the  expense  of  degradations  in 
estimator  performance.  In  passive-sonar  bearing  estimation  problems,  an 
intuitively  logical  approach  is  the  extraction  of  source  bearing  information 
from  the  travel-time  delay  measured  at  the  peak  of  the  cross-correlation 
between  two  spatially  separated  receivers.  Carter  [51  gives  a detailed 
discussion  of  this  sub-optimum  technique  for  a generalized  class  of  cross- 
correlation  processors.  A popular  extension  of  this  approach  to  multiple 
hydrophone  arrays  is  the  split-beam  correlator,  or  tracker  [2],  [18],  [19], 
[36].  This  processor  has  the  intuitive  appeal  of  a cross-correlation  type 
bearing  estimator  and,  under  certain  circumstances,  it  can  rival  the  per- 
formance of  the  optimum  estimator  [19],  For  these  reasons,  we  have  chosen 
the  split-beam  tracker  ns  a representative,  frozen-design,  sub-optimal 
bearing  estimator.  In  later  chapters,  we  analyze  the  detailed  effects  of 
highly  anisotropic  noise  on  the  split-beam  tracker. 

In  section  l of  this  chapter,  we  define  the  split-beam  tracker  (SBT) 
and  briefly  discuss  its  operation.  In  section  2,  we  derive  an  expression 
for  the  SBT  output  in  terms  of  the  array  input  waveforms  and  calculate  its 
mean  and  variance.  In  section  3,  an  approximate  expression  is  developed 
for  the  variance  of  bearing  estimates  generated  by  the  SBT  for  the  particu- 
lar s ignal-noise- Inter ference  model  defined  previously  in  section  2,2. 
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3.1  - The  Split-Beam  Tracker 

A block  diagram  of  the  SBT  array  processor  is  shown  in  figure  3.1. 
The  linear  receiving  array  consists  of  L,  uniformly  spaced  hydrophones 
(mutual  separation  of  d feet).  The  direction  of  maximum  array  sensitivity 
(i.e.  its  steering  angle)  is  controlled  automatically  or  manually  through 
the  selection  of  the  L time-delays,  Tj,  ....  t . The  array  outputs  are  then 
split  in  half  and  summed  to  form  the  left  and  right  half-beam  signals  y^t) 
and  y2(().  The  left  half-beam  output  y^t)  is  phase  shifted  90°  by  time 
differentiation.  Both  signals,  y} (t)  and  y2(t),  are  then  passed  through 
identical  filters  with  transfer  functions  J H ( j oo) | 2 , multiplied  together, 
and  integrated  for  T seconds.  The  processor  output,  z(a),  is  finally 
operated  upon  by  a null  detector  to  obtain  the  target  bearing  estimate,  0. 

The  90°  phase  shift  of  one  split  beam  output  (y^(t))  results  in  a 
mean  processor  output  equal  to  zero  when  the  array  is  steered  exactly  on 
target.  In  the  absence  of  the  90°  phase  shift,  the  correlation,  *(<*),  of 
the  two  split-beam  signals  would  be  a peak  rather  than  a null  when  steered 
at  the  target.  As  the  array  steering  angle,  a,  is  swept  through  180  degrees 
by  adjusting  the  delays  (Tp  t^) , the  mean  processor  output  for  a 

single  signal  in  isotropic  noise  would  increase  from  near-zero  to  a posi- 
tive peak-level  just  prior  to  crossing  zero  at  a = 0^,  followed  by  a negative 
peak  and  a subsequent  return  to  zero-level.  A nominal,  mean  SBT  output 
lor  steering  angles  in  the  region  of  0-180  degrees  is  shown  in  figure  3.2. 

The  task  that  the  estimator  (or  null  detector)  must  perform  to 
obtain  an  individual  bearing  estimate  is  not  quite  as  simple  as  might  be 
inferred  from  figure  3.2.  Since  the  integrator  output,  z(o),  is  a random 

A 

process,  the  source  bearing  estimate,  0,  will  be  a random  variable.  Conse- 
quently, In  practice,  the  algorithm  not  only  estimates  the  angular  location 


ase 


Figure  3.1  - Block  Diagram  of  the  Split-Beam  Tracker 
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of  a particular  zero-crossing,  it  also  decides  between  multiple  zero- 
crossings  and  discriminates  between  positive  and  negative-sloped  zero- 
crossings.  In  theory,  however,  we  can  develop  expressions  for  the  SBT  out- 
put mean  and  variance  and  analyze  its  bearing  estimation  performance  with- 
out actually  confronting  these  problems.  We  demonstrate  this  in  the 
remaining  sections  of  this  chapter. 


3.2  - SBT  Output  Mean  and  Variance 

We  again  assume  that  the  hydrophone  output  signals  Xj(t),  ...»  x^(t) 
are  band-limited,  zero-mean,  Gaussian  processes.  The  outputs  of  the  two 


sp I I t-beamformers  in  figure  3.1  are  simply 

L/ 2 

y, (t ) = £ x (t  + t^) 


1 


e-i 


and 


y^(t)  - y x (t  +t^) 


L 

1 

*=l/2+1 

where  the  steering  delays  t ^ are  given  by 


t ^ = (£  - 1)  • cos  a - (L  - i) 


(3-la) 


(3-lb) 


( ) cos  a 

c 


(3-lc) 


and  where  d Is  the  hydrophone  separation, 

c is  the  speed  of  sound  in  the  sea,  and 

<*  is  the  steering  angle  of  the  array's  main  response  axis. 
Since  individual  SBT  bearing  estimates  are  constructed  solely  on  the  basis 
of  Information  contained  in  a T-second  history  of  y^(t)  and  y^(t),  we  can 
characterize  them  by  their  Fourier  series  representations: 


y.(t)  = (-y-)  1 Y (Jug) 
i=-N 


Jv 


e 


(3-2a) 


36 


,T  1 k 

when-  Yj(j'“n)  “ / yt<t)  e d t i = 1,  2 (3-2b) 

o 

.iml  2'.'n  Miu  total  bandwidth  of  tin-  processes.  Equations  (3-la,  b,  c) 

and  (T-2n,  b)  del  Ine  the  split-beam  output  signals  in  terms  of  the  L hydro- 
phone input-signals. 

The  unsmoothed  processor  output  z(t,a)  is 


z(l,u)  = yj (t)  ’ y ^ ( t ) 

N +ju>Rt\  / 1 


■ (ti’i-.'  •*,v) 


Since  the  transfer  function  of  a time-differentiator  is  H^(jw)  *=  ju>,  the 
unsmoothed  processor  output  can  be  written  as 

N N +j(u>k  + ujt 


x ( t . < 0 -p  £ l (jwk)H(  i<*»k)Yj  (ju)k)H(ja)in)Y2(ju)m)e 


k=-N  in=-N 

The  output  of  the  SBT  integrator  is  simply 

t T 


(3-3) 


(u)  = / z ( t , u)dt 


N N 


k=-N  m=-N  o 


+j((U  + o>  )t 

T k ra  « /o/\ 

e dt  (3-4) 


Now,  if  the  observation  time  T is  much  larger  than  the  correlation  time  of 
the  band-limited  input  processes,  the  integral  in  (3-4)  approximates  the 
Kronecker  delta;  that  is 


/ 


dt  - T6 


k,-m 


(3-5) 


where  T is  the  observation  interval,  and 

| l k = -m 


k , -m 


|o 


k ^ -m 
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Hence , i 1 1 i.-j t f on  (3-4)  and  (3-5)  can  be  combined  to  give,  as  the  SBT  proces- 


sor Ollt  |>I|(, 


1 N 2 

Zo(,l)  = ~ A |H(jouk)|  Yt  (ju>k)Y2(ju.u)  (3-6) 


where  * here  denotes  the  complex  conjugate. 


'•  i om  equation  (3-6)  it  is  an  easy  matter  to  generate  the  mean  SBT 
output  expression: 

V“>  y'J  (jv'H(jwk)'?ElYi(jwk)Y2*(ju)k) I 


T~ J <J“k)|H(Jo,k)|2T  ci2(a,k)  • 

k=-N 


(3-7) 


<:iy(',,k)  is  the  cross-power  spectral  density  between  yj(t)  and  yj(t)  at  the 
i requein  v . Since  the  frequency  summation  in  (3-6)  is  symmetric  about 
N - 0,  we  can  replace  (•^(m^)  by  its  odd  component  (co^  also  being  an  odd 


line  t i on  ) , 


IT  -jw,T  ) 

ODD  fi  (m ,)l=  ODD  )/  R (x)  e dx 


id  i 12k 
k I 


<i».  ) 1? 

k ( o 


fT 

= -j  J K j ( i ) sin  widT 


= j Im  If)  (id  ) ^ 
J | 12  k j 


so  that  ( 1-7)  becomes 


/(,(“)  = ~ 1.  “»k | U(Jn»k>  | Im  jtl12(^k) 

k=-N 

provideil  the  observation  interval  T is  large. 


(3-8) 
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I' rum  (J-j),  we  can  also  obtain  the  mean  square  value  of  the 


SBT 


oil  (put: 


N N 


(3-9) 


IL  is  known  I 18),  that  the  following  identity  exists  for  the  expectation  of 
the  plod  lie t of  complex,  Gaussian  random  variables  required  by  equation  (3—9); 


Y 

I I 


( i o,k ) Y*  ( j oik ) Y* ( j fom) Y2  ( j «om) j = E|Y1(jWk)Y*(ji,k)|  E j Y*^^) Y2 (Jui^) j 

+ E|Y,  (ja)k)Y*(jum)j  EjY^J^Y^Jw^j 

(3-10) 

I lie  I Irsl  term  c> f (3—10)  is  simply  the  product  of  the  cross— power  spectral 
density  at  two  different  Irequencies,  i.e. 

EjV'VVH)!  E !Y*<J“m)Y2(J“»5J  = t2gi2H)G12*(ub>  (a-U*) 

Relative  to  the  second  terra  of  (3-10),  consider  the  following: 

fT 


l':Jvi(i"k>vl(j"„  | ■ r dt  /'  du  F Jy,  (t)y,  (u)j  e 3<“k‘  “”U> 


./* -....to  ;lvf;,(v->' 


du 


= T Cll(l‘'k)6k,m 


r 


T T I \ ~j(u>,  t+  <*>  u) 

K Yl(J«k)Y2(J“«)|  = / dt  / du  E y,(t)y2(u)'  e 


o o 


= j <JtR12(t)  e / 


JW.T  ~j(u)  + U>  )ll 


du 


= T C (w,  ) 6. 

12  k'  k,-ra 

Consequent  ly,  tl»e  last  term  of  (3-10)  can  be  written  as, 

F.jv,(J»k)V!(J„|<)j  ' T2G,2(^)ct2(Um)«ki.m 

(3-llc) 

t 

so  lli.u,  upon  substitution  of  (3-10)  and  (3-lla,  b,  c)  into  (3-9),  the  mean 
square  SBT  outp  t becomes, 

N N 

5 

k= 


l.l<  /.  I. 

I o o , 


l l )|H(>  )|!|ho»  )|2jr2c12<^)r'2(u  ) 

=-N  m*-N  ( 

+ T G,  , (alk)G22(uik)6k^m  + T Gj  2 '<VGl  2 (u)m>  6R,-m  | 


= LUJulk)'H<H)l  Gi2(a>k)]  '[j  <J“B>IH(J«B>I  G12(V]* 

?2  4 N 2 4 * 

uiJH(*u’k*  I G1 1 ^a,jc^G22fwic)  + 1 | HC  | Gj  2 (Wjc)Gl  2 ^ 

k=*-N  k=-N 

(3-12) 

From  oqu.it  ions  (3-7)  and  (3-8)  It  is  obvious  that  the  first  term  of  (3-12) 

— — — 

Is  equivalent  to  z • z = z • z.  Therefore,  we  can  write  the  following  ex- 


pression lor  the  variance  of  SBT  output; 
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V"’  ■tLN“k|,1<l^)l  ■ -'wv  !•  «-i3) 

MacDonald  and  Schul theiss  [19]  have  previously  developed  similar 
expressions  for  the  SBT  output  mean  and  variance  given  by  equations  (3-8) 

and  < 3— 13).  Discussion  of  the  form  of  the  split-beam  filter  transfer  func- 

2 


tious,  !ll(Ju>k)|,  is  deferred  until  the  end  of  this  chapter.  At  this  point, 
vc  can  apply  the  general  expressions  in  (3-8)  and  (3-13)  to  the  specific 
;lgnal-noise- interference  model  of  interest.  Recall  from  chapter  2 that 
i lie  output  of  the  Cth  hydrophone  can  be  written  as 


M 

xf(t)  = s(t  - A£)  + n£(t)  + l im(t  - 6®)  (3-14) 


i lie  required  spectral  densities  are,  by  definition 


T 

71  -JV 

Gi  i = Rii  ^ e dT 

J -T  3 


i,  j = 1,2 


(3-15) 


where  K . . ( i ) is  the  correlation  function  of  y^(t)  and  (t  + t)  for  i,  j * 1,2 
and  where,  without  loss  of  generality,  we  redefined  the  observation  interval 


The  correlation  functions  in  (3-15)  are,  by  definition. 


RjjU  - u)  = h J y ± ( t ) yj  (u)  J 


(3-16) 


from  tin's  definition  and  equations  (3-la)  and  (3-14)  the  autocorrelation  of 
,(t)  is, 
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R,  1 1 0 = e j yt  (t)yj  <u) 


li  I | 

e|(L  “«(t  4 'f’pL  vu  4 vj 


l l eJ" js(t  - A.  + T„)S(u  - A + t ) J + Ejn«(t)n  (u)j 

J-l  p-1  U 4 4 p p > * 4 I 

l E | i (t  - S"  + x„)i  (u  - 6m  + t ) J 1 
I m 4 1 m P P I J 


(3-17) 


where  we  have  made  use  of  the  fact  that  the  signal,  noise  and  all  inter- 
ference processes  are  mutually  independent.  We  can  further  simplify  (3-17) 


by  writ i ng 


L L 
2 2 


R,,(i)  = l l KQB(T  - a;  + A')  + R (t)6 » 
1*1  p=l ' 4 P nn  cp 


* I Rm.(t  - a + v )l 

m=l  11  Zn  pm  I 


where  = A„  - t 0 * — (cos  9_  - cos  « ) + — cos  a 

■c  4 4 c T c 


(3-18) 


jl  l = p 

fp  |o  l # p 


, sl. -sn  rcos „ . cos  „ 1 +|-<k^w  cos  a I 

tm  c [ ra  J L c j 

Substitution  of  (3-18)  into  (3-15)  yields,  after  some  algebra, 
L L 

2 2 ( -Jw  (Ai  - A') 

Gl'(,‘,k)  m l \ ) S(“k)  e + N(\)6lp 

4=1  p=l  I r 


+ t Im(wk)  6 

m=l 


-n<f’in  - S-> 


(3-19) 
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wiiere  S • *>k) » NG^),  ond  I (u>  > are  the  signal,  noise  and  mth  Interference 
power  spectral  densities,  respectively,  at  the  frequency  In  a similar 

fashion  it  Is  possible  to  obtain 

L 


L L ( -jw  (Al  - A') 

l l )S<V  e k 1 P + N(u,k)6£p 


C?2(u’k)  | f'“k 

1 P^+i 


(3-20) 


M -j w.(Ci  - V )\ 

+ j , ( , . * p« 

m=l 

The  cross-correlation  between  (t)  and  y2(t)  is 

r12(t)  = E}y1(t)y2(u)  | 

L 

2 L r 

= l L Pjs(t  - A,  + T «) s (u  - A + t )}  + E n»(t)n  (u)| 
1*1  p=^-lL  < L P P j I L P » 

M 


\LEi1"<t~{"  + l£>i'»<“'5"  + Vl] 


L 

2 L 


■ I L ks<T  - a;  + a:>  + 1 O*  - cl,  + e* )] 

e-i  p=4+l(  33  ^ p m=l  11  Pm  j 


Consequently,  the  cross-power  spectral  density  G12(wk)  is 

L 

2 L ( -joa  (Ai  - A’)  M 

G.2(“k>  -J,  M?<V  e p + 1 W 


£=1  p=»2+l 


- V] 


(3-21) 

where  A£,  Ap,  C^,  and  £'m  are  defined  as  in  equation  (3-18). 

Substituting  equation  (3-21)  into  equation  (3-8)  yields  the  following 
express! "ii  of  the  mean  SBT  output  for  the  signal-noise- interference  model 


of  interi  st : 
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N 2 f 2 L 

zo(a)  =,  K,  (i,k|H(jwk)i  l k 


k=-N 


l- 1 P«f  +1 


sK)[sin  “k(Ai " ^3 


+Jl  Sin  [Wk(4m  " ^ 


(3-22) 


Likewise.  the  substitution  of  equations  (3-19),  (3-20)  and  (3-21)  into 
(3-13)  yields  the  following  expression  for  the  SBT  output  variances: 


2 N 2 4 

’ ('«)  = I “I  |H(j(.)  ) | 

Zo  k=-N  k 


M -Jo).  (€'  " 4’  ) 

3-  l I U.  ) e k Pm  qm 

u m k 


m=l 


-jw  (A'  - A') 

S(u.  ) e k r s + N(4).  )6 
k k rs 


+ ? . f , -n(^i  - «ii> 

+ ),  1.(01  ) e 

i=l  1 K 


M -jw.U'  " V )'  2 

+ l l (a,,)  e k Pm  qm 
4 m k 


in=l 


L 

2 L -jut,  (A'  - A’) 

k P q 

L,,~'~k'  - 

T 


(3-23) 


l lh  SK>  e 
lp=l  q=^r+l 


7 L L 

\ 2 2 L L 

-jm  (A*  - A') 

hill 

1 L L 

S (w.  ) e R P q + N(u.  ) 6 

k k pq 

^>=1  q=l  r=*2+l  s=yH 

i 

3.3  - Theoretical  SBT  Bearing  Estimation  Variance 

figure  3.3  shows  that  the  steering-angle  (a)  dependence  of  the 
theoretical  mean  SBT  output  (z^(a)  in  equation  (3-22))  closely  follows  the 
eon jectured  u-dependenee  of  the  nominal  SBT  output  in  figure  3.2.  The 
effect  ol  a single,  weak  interference  at  45°  is  evident  in  figure  3.3  while 
the  effect  of  the  target,  located  at  90°,  is  most  obvious.  The  linearity 
of  («)  over  the  small  angular  region  about  0^,  the  true  source  bearing, 
provides  the  basis  for  a theoretical  approximation  to  the  bearing  variance, 
11  , of  the  SBT  processor.  MacDonald  and  Schultheiss  [19]  first  reported 


2 _ 

a straightforward  approximation  to  o.  that  assumes  z (a)  is  linear  over 
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at  1 eusL  an  <«  = 2o^  region  about  0T-  The  nature  of  their  approximation  can 

he  clearly  seen  in  figure  3.4.  If  z^fct)  is  linear  for  “ in  the  region 

|0  ~ n , 0 + o | and  if  a (u)  is  relatively  insensitive  to  changes  in  a 
o o z 

o 

over  this  region,  then  the  following  linear,  stochastic  model  of  z (a)  is 
vii  lid: 


z (“)  = A + Ba  + £,  G-o  < a < 0 + a 

o e q 


(3-24) 


where  t.  Is  a zero-mean,  Gaussian  random  variable  with  variance  equal  to 

<i  . The  estimate  G occurs  at  the  zero-crossing  of  z (a)  so  that, 

'o  0 

- 0 - a + ■» + < . 

Krom  tills  model,  we  obtain 


so  that. 


and 


B 


Since  I!  =■  , we  obtain  the  following  expression  for  the  MacDonald- 

Schul thelss  theoretical  approximation  to  the  bearing  variance  of  the  SBT 
processor : 
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o"  (.*) 
z 

o 


(3-25) 


is  unbiased,  then  0 - 0 and  (3-25)  becomes 


(3-26) 


I i Ls  bliised,  but  only  slightly,  then  equation  (3-26)  is  still  a good 

2 

approx imnt Ion  to  o . Otherwise,  equation  (3-25)  is  more  accurate.  If, 

_ 2 

however,  one  can  assume  that  z («)  is  linear  and  a (a)  is  almost  constant 

o z 

. o 

o i r a larjjejr  angular  region  about  G,  then  equation  (3-26)  is  again  valid. 

I’,  for  example,  the  true  target  were  located  at  the  bearing  in  figure  3.4, 
i 'ion  (3-26)  is  equally  as  valid  as  (3-25)  since  the  assumptions  are  true  for 


" (»  u ,,  • 

The  error  metric  defined  in  (3-26),  although  not  as  general  as  that 
ii  (3-25),  Is  somewhat  simpler  to  interpret  since  it  decouples  the  variance 
aid  bias  errors.  Equation  (3-26)  is  the  approximate  theoretical  SBT  variance 
in,  trie  examined  in  chapters  5 and  6.  The  numerator  and  denominator  of 
e,  nation  (3-26)  exist  in  the  form  of  equations  (3-22)  and  (3-23).  Equation 
( -.12)  must  be  differentiated  with  respect  to  a and  both  equations  then 
o'.ilunted  at  = G^,.  Recalling  the  definitions  of  Aj,and  , the  derivative 
o'  I lie  mean  SBT  output  can  be  written  as. 


_k 
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3 / N 
* 1 


l vl»Ou>| 

;i"  k— N 


L 

2 I. 

X L 

H P-  f+l 


(sin^(wk)  cos  [wk(Aj,  - A’)] 


' i wi 

tn  - 1 l< 


-- sin“  1 (u),  ) cos  fu>.  (4)  - V )1 

c m k L k In  pm  J 


(3-27) 


To  ova  I n ite  (3-27)  at  a = 13  , note  that 


A'  = 0 
l» 


Hill 


'•In 


- r ’ = — [cos  0 - cos  0 ] 

pn  c mi 


(3-28a) 


(3-28b) 


so  that  aquation  (3-27)  becomes, 


i)z  (“) 
o 


d sin  0 N 2 ?(  ^ L r- 

= I w |H(jw  )|  1 l fL  (£  - p)S(u>k) 

c k=-N  k Vf-l  P=S+1L 

U r 0 2 


T 


M 


+ U - p)  l xm(wk)  cos  [(£  - p) si] 


m=l 


i 


(3-29) 


where  we  have  defined  U = u.  -—[cos  0 - cos  0 ].  As  shown  in  appendix  C, 

k c m l 

|.v  expressing  the  cosine  term  in  (3-29)  as  a sum  of  complex  exponentials, 

(i  (s  possible  to  evaluate  both  double  summations  over  the  hydrophone  indices 
(i.e.  equations  (C-J)  and  (C-8))  and  obtain  the  following  expression  for 
the  denominator  of  (3-26) 


0z  ( •*) 
o 


u = 


c , k1 
k=-N 


k>i2| 


-d  sin  0T  N 2 J M 


l -klH<>k>l  li  s(“k>  - 1 ww 

ra=l 


x £sin(L  ^)  cot(-&)  - j sin(L^)  cot(|  • |)  - | cos  (L|)jj 


(3-30) 
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where 


C (w.  ) 
m k 


2 i tt 

•sin  q ~) 


•4  in  (—) 


and 


■i  = m. 


d r 

cos 

c L 


0 - cos  6 

m i 


Using  the  results  of  appendix  C and  equations  (2-28a,  b) , it  can  be  shov;n 
that 


o CO 

z 


t=e„ 


N ^ it  ( 2 ? J 

= l u)k|H(jo)k>|  <-£  N (u)k)  + ^ S(uk)N(wk) 


k=-N 


M 


+ LN<“k>  1,  WW  + L"SK>  " VVW51"^  • 

m=l 


+ 4 


r M 

\i 

lm=l 


M L n - r M 

I (u>.  )C  (to,  )sin£  ~) 
m k m k 2 2 


-i  r M 

•Mi  (u  )c 

J L =i  n k n 


(a)  )cos(  j 


k' 2 

(3-31) 

The  final  expression  for  the  approximate,  theoretical  SBT  bearing  variance 
is  given  by  equation  (3-26)  with  its  numerator  and  denominator  defined  ex- 
plicitly by  equations  (3-31)  and  (3-30),  respectively.  This  performance 
metric  will  be  evaluated  numerically  in  chapters  5 and  6. 

I'he  bearing  variance  metric  presented  above  is  expressed  in  terms 

2 

of  the  filter  transfer  functions,  |H(jwk) | , which  are  as  yet  unspecified. 
Ideally,  one  should  like  to  select  that  filter  transfer  function  which  min- 
imizes the  bearing  estimation  variance.  Since  the  sub-optimal  SBT  processor 
is  designed  to  discriminate  only  against  isotropic,  Gaussian  noise,  j H ( j ) I 
cannot  bo  selected  by  minimizing  the  above  bearing  variance  expressions 
(equations  (3-26),  (3-30),  and  (3-31))  which  pertain  to  a highly  anisotropic 


* 


noise  model.  However,  if  we  let  I (u>.  ) = 0 for  m = 1,  ....  M at  all 

m k 

frequencies,  then  our  bearing  variance  expressions  will  reflect  a noise 
model  which  is  isotropic.  Under  this  isotropic  assumption  (3-31)  and  (3-30) 
can  be  substituted  into  (3-26)  to  obtain: 


It  is  known  fl8],  [19]  that  the  bearing  variance  in  (3-32)  can  be  minimized 
through  a variational  calculus  approach  and  results  in  the  following  filter 


funct ion: 


|HU“’k) 


S(mk) 

N*('uik) 


1 + 


LS(mk) 

N(mk) 


(3-33) 


Use  of  this  filter  transfer  function  in  the  SBT  processor  shown  in  figure 
3.1  will  optimize  its  bearing  estimation  performance  for  a unidirectional, 
plane  wave,  Gaussian  signal  in  isotropic,  Gaussian  noise.  In  the  remainder 
of  this  research,  when  a specific  form  of  the  SBT  filter  coefficients  is 
required,  the  filter  defined  by  (3-33)  will  be  employed. 


CHAPTER  IV 


DIGITAL  SBT  SIMULATION 

Tills  chapter  describes  a digital  computer  simulation  which  was 
developed  to  provide  additional  measures  of  the  performance  of  the  SBT  in 
highly  anisotropic  noise  environments.  A real-time  array  processor  can  be 
Implemented  in  either  the  time  or  frequency  domains.  Traditionally,  sonar 
processors  have  been  implemented  in  the  temporal  domain.  However,  with  the 
current,  sophisticated  digital  technology,  sonar  processors  are  also  easily 
Implemented  in  the  frequency  domain. 

Discrete  models  of  some  temporal  operations,  such  as  differentia- 
tion, ire  unavoidably  noisy  procedures.  However,  the  equivalent  frequency 
domain  operation  is  complex  multiplication  and  easily  modeled.  There 
are  also  problems  associated  with  digital  implementations  or  simulation 
models  of  array  beam  steering.  Discrete  temporal  beam  steering  requires 
extremely  high  sampling  rates  or  the  use  of  complex  interpolation  algorithms. 
In  the  frequency  domain,  beam  steering  is  also  a simpler  operation.  At 
each  frequency,  the  input  is  simply  multiplied  by  a complex-exponential 
weighting.  The  advantages  of  frequency  domain  implementations  (or  simula- 
tions) are  often  completely  negated  by  the  added  cost  and  complexity  of 
the  t iine-f requency  transformation.  If  the  input  process  spectra  vary 
considirably  with  frequency  and  exist  over  a large  bandwidth,  then  even 
the  computational  speed  of  FFT  algorithms  may  fall  short  of  making  fre- 
quency domain  processing  attractive. 


In  section  4.1  below,  we  discuss  the  reasons  for  developing  a 
frequency  domain  simulation  and  describe  its  operation.  Section  4.2 
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describes  the  methodology  employed  to  generate  samples  of  the  desired  in- 
put, slgnal-nolse-interference  processes  at  each  hydrophone.  The  last 
section  develops  the  digital  bearing  estimation  algorithm. 


4.1  - Digital  SBT  Simulation 

A block  diagram  of  the  complete  frequency-domain  SBT  simulation 
is  shown  in  figure  4.1.  The  simulation  performs  three  basic  functions: 
generation  and  FFT  processing  of  the  input  data  sequences,  application 
of  the  SBT  algorithm,  and  extraction  of  the  source  bearing  estimate  from 
the  processor  output  The  SBT  as  configured  in  figure  4.1  is  conceptually 
capable  of  processing  recorded  field  data  instead  of  the  T-second  long 
simulated  time  sequences,  x^(nAt),  £* 1,  ...,  L.  However,  for  the  purpose 
of  emulating  the  SBT  processor,  observe  that  if  the  frequency  domain 
characterization  of  the  hydrophone  outputs  could  be  synthesized  directly, 
then  the  costly  requirement  for  a Fourier  transformer  is  eliminated. 
Indeed,  the  linearity  of  the  Fourier  transform  assures  that  the  Fourier 
coefficients  wi LI  be  Gaussian  random  variables  when  the  hydrophone  outputs 
are  Gaussian  processes.  Consequently,  the  direct  generation  of  the 
frequency-domain  random  deviates  is  justified.  In  the  next  section,  we 
describe  a methodology  for  generating  the  Gaussian  deviates, 
k-0,  1,  ...,  N-l,  which  possess  the  desired  statistical  properties. 

The  hydrophone  output  sequences  of  random  deviates  stimulate  the 
SBT  array  processing  algorithm  which  generates  an  output  random  deviate 
corresponding  to  the  chosen  steering  angle,  a.  If  the  steering  angle  is 
subsequently  incremented  and  the  SBT  algorithm  reapplied  a total  of  ISTEER 
times,  then  the  processor  output  sequence  z^O),  zq(2),  ...,  zq(ISTEER)  is 
generated.  The  simulated  bearing  estimates  are  then  formed  from  these 
ISTEER  deviates. 


V U <J 

c a c 


u 

Q> 

u 

u 

oc 

cn 

£ 

3 u 

c 

rH 

a» 

Cu  o 

•H 

<0 

3 

B ** 

CO 

£ 

O' 

o 

(0  4J 

*0 

0> 

U E 

D £ 

x: 

U O 

w 

U-i 

T~\  U 

0 X 

CX3  *H 

*-■  a 

4J 

0/  £ 

4-*  U 

fX  o 

3 

x:  »h 

•H  O 

>-« 

a. 

VJ  £ 

QC  aoh  T3 

4-1 

B 

•H  r— 4 

PQ  >. 

3 

£ O 

c < 

C/)  ,£ 

O 

•H  T3 

1 

£ 

0) 

1 

•H 

i) 

Vi 

• 

4-J 

TJ 

U-i 

V 

-A 

0) 

CN 

(0 

C 

o 

u 

(ft 

•H 

<0 

£ 

XJ 

*H 

U-l 

> 

£ 

•H 

£ 

o 

O 

y 

£ 

o 

£ 

TJ 

O 

•H 

2: 

if) 

•H 

4J 

• 

4-1 

0) 

£ 

4J 

<0 

£ 

u 

£ 

(0 

£ 

CJ 

y 

£ 

•H 

u 

•H 

JO 

£ 

u 

0) 

to 

V 

£ 

a) 

V 

3 

CO 

si 

£ 

e 

OC 

v 

cr 

3 

u 

0) 

0 

•H 

a> 

0 

a 

JO 

(0 

o 

u 

C/5 

Uh 

u 

co 

o 

a> 

Frequency  Domain  SBT  Simulation 
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The  SKT  processing  algorithm  is  simply  the  output  expression 

developed  in  . h.ipter  3 tor  arbitrary  hydrophone  output  waveforms  and  need 

not  be  discussed  in  detail  again.  Specifically,  the  SBT  algorithm  is 

defined  by  equations  (3-1)  through  (3-3)  and  equation  (3-6)  where  the 
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filter  transfer  function  |H(Jcd^)!  is  defined  in  equation  (3-33).  Recall 
that  the  utility  of  the  frequency  domain  processor-simulation  arises  from 
the  ease  of  performing  the  following  operations  in  the  frequency  domain: 
beamsteering,  differentiation  (or  90°  relative  phase  shift),  filtering  and 
Integrat ion. 

The  third  function  performed  by  the  simulation  is  the  actual  bear- 
ing estimation  algorithm.  We  utilize  a linear,  null-detection  algorithm 
which  closely  follows  the  one  developed  by  Witt  [36],  That  is,  the  bearing 
estimate  is  the  angular  location  of  the  steepest,  negatively-sloped,  zero- 
crossing  of  the  linear  regression  of  the  output  sequence,  zo(l),  ..., 
zq(ISTKER),  on  the  steering  angle  sequence  tXj,  ....  ajcpEER"  0ur  a^8orir-hm 
is  detailed  in  section  4.3. 

4.2-  Generation  of  Input  Random  Deviates 

Before  discussing  the  direct  generation  of  the  frequency  deviates, 
X^CJoJr).  it  is  of  interest  to  outline  the  steps  required  to  generate  the 
tLme  sequences,  x^(nAt).  Since  the  noise  model  of  interest  contains  M- 
directional  interferences  and  an  uncorrelated,  isotropic  component,  a total 
of  M + L + 1 independent  sequences  of  Gaussian  random  deviates  are  required. 
At  each  hydrophone,  the  M+l  sequences  associated  with  the  target  and  Inter- 
ferences must  be  time  delayed  in  proportion  to  their  respective  bearings. 

As  mentioned  earlier,  this  time  delay  insertion  requires  high  sample  rates 
or  interpolation  between  every  adjacent  point  in  each  of  the  L hydrophone 
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l line  sequences.  The  M + I-  + 1 sequences  are  discrete  samples  of  white, 
Gaussian  noise  processes  and  must  each  be  filtered  to  obtain  the  desired 
spectral  characteristics.  If,  for  example,  one  restricts  his  attention  to 
low-pass  band- 1 i mi  ted , white  Gaussian  processes,  then  the  spectral 
shaping  is  easily  accomplished  by  a linear-phase,  FIR  digital  filter  [24]. 
Finally,  the  simulated  output  signal  at  each  hydrophone  is  the  sum  of  an 
independent  noise  sequence  and  M + 1 properly  delayed,  propagating 
(directional)  sequences.  As  indicated  in  figure  4.1,  the  input  time 
sequences  must  also  be  pre-processed  by  the  Fourier  transformer  in  order  to 
stimulate  a frequency-domain  SBT  algorithm.  If  one  intends  to  stimulate 
a time  domain  algorithm,  it  has  been  noted  [36]  that  the  input  time  se- 
quences must  be  of  length  greater  than  T to  permit  truncation  of  noisy  end- 
point values  generated  by  the  time-differentiation  operation. 

If,  as  is  the  case  here,  one  can  define  the  first  and  second  order 
signal-plus-total-noise  statistics,  then  it  is  possible  to  generate  the  re- 
quired Input  random  deviates  through  a transformation  of  zero-mean,  unit 
variance  random  deviates  (e.g.  [4]  and  [23]).  Our  approach  closely 
parallels  that  of  Franklin  [10].  By  assumption,  the  hydrophone  outputs 
are  zero-mean,  Gaussian  random  processes.  By  design,  the  hydrophone  output 
Xf(jV  at  frequency  and  the  output  ^(joi  ) at  frequency  u>m  are  independ- 
ent for  iuk  # i»m.  Consequently,  the  array  response  is  completely  specified 
by  the  KK,  L x L dimensional  signal-plus-total-noise  covariance  matrices, 

K„„.  In  other  words,  for  k = 0,  1,  ...,  KK-1,  we  have  at  the  frequency  u> 

— oN  K 
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X = j X t ( j !.>k)  , X^(ja)k),  ....  XjC.j  0)k>  1 


»•:  i x o 


E'  XX  } = KCHU.  ) 
( I — SN  k 


The  covariance  matrix  Kc  (w  ) is  not  only  Hermitian,  it  is  also  Toeplitz 

bN  K 

and  Is  defined  In  equation  (3-20)  as. 


■W  * TS<"k>  i("k>  • *<“k>  + ™ <Vi+Il  W V“k>  C(V 

m=l 


-X  * 

— TS  a a + 


(4-D 


where  S,  N,  and  1 are,  respectively,  the  signal,  noise,  and  interference 

-A  -A 

l»ower  spectral  densities  at  uo^,  a (w^)  and  b^  (w^)  are  the  signal  and  mth 

Interference  phase  delay  vectors,  T is  the  observation  time  interval,  and 

K„  is  the  noise  plus  interference  covariance  matrix  at  o>,  . Since  both  the 
--N  k 

propagating  and  non-propagating  components  of  the  received  waveforms  are 
by  assumption  homogeneous,  the  signal-plus-total-noise  covariance  matrix 
in  (4-1)  is  a Toeplitz  matrix.  For  simplicity,  we  drop  the  frequency  index 
for  the  remainder  of  the  derivation. 

-A 

Generation  of  the  desired  complex  vector  Gaussian  deviate  X first 
requires  the  generation  of  an  L dimensional,  complex  vector  Gaussian  variable 
7 which  has  zero  mean  and  unit  variance.  The  vector  z is  easily  constructed 
from  two  real,  Gaussian,  vector  random  variables  with  zero  mean  and  unit 


variance  . If  the  probability  density  functions  of  the  real,  independent, 
l.-d  Imens ional  vectors  u and  v are  N(0,  _I) , then  by  defining 


2',  (i)  [?+Jvj 


(4-'») 
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it  is  apparent  that 


I I - 
i:  i as  i - l) 


K z z 


A standard  DKG  system  KI.-10  library  subroutine  was  used  to  generate  the 


zero-mean,  unit  variance  Gaussian  elements  of  the  vector  random  variables 


u and  v. 


In  general,  the  vector  X is  related  to  the  vector  z by  the  follow- 


ing transformation: 


X = C z + u 


(4-3) 


where  C is  a unique,  L x L lower-triangular,  deterministic,  complex  matrix 


such  that 


Ksn  = E I x x [ = c • c 


(4-4) 


and,  by  assumption,  p = E j X J = 0.  Since  is  known  (see  equation  (4-1)), 


equation  (4-4)  can  be  solved  recursively  (Franklin  [10],  p.  208)  for  the 


elements  of  the  transformation  matrix  C. 


The  signal,  noise,  and  interference  power  spectral  densities  must 


be  defined  at  each  frequency  in  order  to  completely  specify  K . To 


facilitate  numerical  evaluation  of  both  the  theoretical  and  simulation 


error  measures,  we  make  the  simplifying  assumption  that  the  spectral  shapes 


oi  all  processes  are  identical,  i.e. 


SU),  N(u),  I (u»)  = (S,  N,  1 ) • G (u.) 
m in 


(4-5) 
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where  , N anil  1 1(  in  (4-5)  are  constants  and  C(m)  is  a normalized  power 
spectr  I dens  i i function.  In  the  remainder  of  this  study,  we  shall 
assume  in  observation  interval  of  T = 0.128  seconds  and  a sampling  fre- 
quency of  250  II/..  Consequently , the  one-sided  spectrum  G ( to)  is  defined 

2 Jl 

by  16  amples  at  k,  k = 0,  1,  2,  , 15.  Furthermore,  we 

requir<  the  low  pass  spectral  sequence,  GO^)  to  roll-off  at  40  dB/octave 
and  del  ine  its  12  dB  down  frequency  to  be  uij  j = ^ • 13  = 211(101.5  Hz). 

Figure  't.2  Illustrates  the  idealistic  spectral  shape  of  signal,  noise  and 
interl.  rence  employed  in  this  study.  Since  the  parameters  — 9^, 

M,  and  ■)  (m  = l,  . . . , M)  are  investigated  as  variables  within  this  study, 
we  havi  completely  specified  Kg^  and  are  able  to  use  the  above  methodology 
(equations  (4-1),  (4-4),  (4-3)  and  (4-2))  to  transform  independent,  zero- 
mean,  unit-variance  Gaussian  deviates  into  the  desired,  complex,  vector 
Gaussian  deviates  X(joj^)  (k  = 0,  1,  ...,  15)  from  the  distribution 
N(0,  KM(u,k)). 

5.3  - S i mu lated  Kstimation  Algorithm 

If,  as  In  chapter  3,  we  make  the  assumption  that  the  mean  SBT  out- 
put is  i linear  function  of  steering  angle,  (around  6),  then  it  seems 
reasonable  to  form  an  LSE  straight  line  fit  to  several  output  points  which 
encompass  a zero  crossing  and  calculate  the  bearing  estimate  from  this  line. 
The  simulation  program  calculates  "ISTEER"  SBT  output  deviates  centered 
around  i he  true  target  bearing,  from  which  the  line  fit  is  computed.  In 
the  baseline  numerical  evaluation  that  follows,  the  parameter  ISTEER  is 
equal  t.*  13.  For  simplicity,  let  P =>  ISTEER  in  the  following  derivation. 


saassss 


Figure  4.3  illustrates  the  estimation  procedure  for  V hypothetical 
simulation  output  values.  The  LSE  line  is  defined  by 


ZI.SE(a)  = A “ Ba 


(4-6) 


The  bearing  estimate,  0,  satisfies 


We>  ■ 0 


and  is  clearly, 


(4-7) 


The  coefficients  A and  B are  chosen  to  minimize  the  mean-square-error. 


P th  (Z±  ’ ZLSE;(ai)) 

■fJi(zi  - (A  “ ’V) 


(4-8) 


By  setting  the  derivatives,  with  respect  to  A and  B,  of  (4-8)  equal  to 


zero,  A and  B are  easily  found  and,  upon  substitution  into  (4-7)  yield, 


l zj  l aj  -l  21ai  l aj 


(4-9) 


where  all  summations  range  from  1 to  F.  Since  we  are  interested  in  target 
tracking  applications,  as  long  as  the  sector  over  which  the  beam  is  swept 
is  not  unrealistically  small,  we  need  not  be  unduly  concerned  with  multiple 
zero-crossings  of  the  simulation  output.  Consequently,  the  simulation 
program  uses  all  of  the  ISTEER  output  values  to  calculate  0 and  requires 
only  that  the  slope  of  z (a)  be  negative. 


Figure  4.3  - Discrete  Bearing  Estimation  by  LSE  Line  Fit 
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Previous  studies  [e.g.,  36]  have  successfully  employed  similar 
techniques  for  estimating  source  bearing  from  sampled  SBT  output  values. 

As  discussed  in  chapter  5,  there  are  potential  problems  with  this  type  of 
algorithm.  It  is  desirable  to  have  a large  number  of  beam  positions, 
ISTEER,  as  well  as  having  the  angular  beam-sweep  region  large  relative  to 
the  hearing  standard  deviation.  Otherwise,  it  is  likely  that  LSE  lines 
with  very  small  slopes  will  occur,  resulting  in  unplausible  bearing 
estimates.  In  any  event,  this  problem  requires  larger  simulation  sample 
sizes. 

For  every  parametric  combination  (SNR,  M,  L,  6 , etc.)  of  interest, 

m 

the  simulation  must  be  exercised  a total  of  ISTEER  x times,  where 

KcJr 

J__D  is  the  number  of  estimates  per  sample  and,  recall,  ISTEER  is  the 

Kfcr 

number  of  beam  positions  per  estimate.  Based  on  these  samples,  statistical 
measures  of  the  performance  of  the  SBT  array  processor  are  defined  in  the 
next  chapter. 


CHAPTER  V 


NUMERICAL  RESULTS  - BASELINE 

This  chapter  presents  the  iesults  of  the  numerical  analysis  cf 
bearing  estimation  in  the  presence  of  multiple,  interfering  noise  sources. 
The  Cramer-Rao  (CR)  bound  derived  in  chapter  2 is  the  performance  metric 
for  examining  the  effects  of  such  interference  on  the  optimal  (ML)  esti- 


i 

mator,  which  is  assumed  to  be  unbiased.  For  the  sub-optimal  split-beam 


I 


tracker  (SBT),  bias  as  well  as  standard  deviation  of  error  is  evaluated. 

The  performance  of  the  SBT  is  presented  in  terms  of  both  theoretical 
measures  of  error  and  statistical  measures  of  the  simulated  SBT  perform- 
ance. The  SBT  performance  metrics  are  described  in  detail  in  the  follow- 
ing section.  Section  2 lists  the  baseline  values  of  the  fixed  and  variable 
array,  processor,  signal,  noise  and  interference  parameters.  Section  3 
presents  and  discusses  the  baseline  results  in  terms  of  the  standard 
deviation-type  error  measures  while  the  last  section  contains  the  baseline 
results  In  terms  of  the  bias  error  metrics. 


5.1  - SBT  Performance  Metrics 

One  theoretical  measure  of  the  SBT  performance  is  the  approximate 
expression  of  the  standard  deviation  of  the  bearing  estimate  derived  in 
chapter  3.  The  approximation  arose  from  the  assumptions  of  long  obser- 
vation Interval  and  linear  behavior  (as  a function  of  steering  angle)  of 
the  SBT  output.  Another  theoretical  measure  of  the  SBT  performance  is  the 
estimation  bias  or  expected  value  of  the  bearing  estimation  error.  It  is 
reasonable  to  assume  that  the  SBT  bias  error,  unlike  that  of  the  maximum 
likelihood  estimator,  will  be  significant  for  many  target-to- interference 
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{operations.  The  bias  error,  B,  is  defined  as 

0 - 


B - e !e | - en 


(5-1) 


where 


e‘* 


I I 


denotes  the  statistical  expectation,  6 denotes  the  bearing 


estimate  and  0T  is  the  true  value  of  the  target  bearing.  Assuming  the 
probability  density  function  of  the  SBT  output,  zq,  is  symmetric  about  its 
mean,  the  expected  value  of  the  bearing  estimate  can  be  calculated  from: 


Eizo(a:!  I “ “ Zo(a)i  " “ 0 (5-2) 

' ' I a=0  I a=0 

where  z («)  is  the  mean  SBT  output  derived  in  chapter  3.  The  exact  solu- 
tion of  equation  (5-2)  for  0 is  mathematically  intractable.  A good 
approximation  of  0 is  obtained  for  forming  a least-square-error  (LSE) 
straight  Line  fit  to  p values  of  ZQ(o)  and  solving  the  result  for  0 accord- 
ing to  equation  (5-2).  This  algorithm  is  the  same  one  employed  to  obtain 
bearing  estimates  from  the  simulated  SBT  output  which  was  described  in 
chapter  4.  The  bias  error  in  equation  (5-1)  above  and  the  standard  devia- 
t Lon  of  the  bearing  estimate  given  by  equation  (3-26)  are  the  theoretical 
measures  of  SBT  performance  which  are  presented  in  this  chapter.  They  both 
represent  approximations  to  ensemble  averages  of  bias  error  and  estimation 
standard  deviation  assuming  a linearity  of  the  SBT  output  with  respect  to 
steering  angle,  a. 

Ihe  relevant  statistical  measures  of  the  performance  of  the  simu- 
lated SBT  processor  are  simply  the  sample  mean  and  standard  deviation  of 
the  estimated  bearings.  Samples  of  sizes  from  20  to  60  replications  were 
generated  for  the  various  parametric  combinations.  Consistent  with  the 
previous  theoretical  measure,  the  sample  bias  error  is  presented  in  this 


■ 
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chapter  rather  than  the  sample  mean  of  the  estimated  bearing.  With  the 

* 4 

1th  bearing  estimate,  0^,  generated  by  the  simulation  as  described  in 
chapter  4,  the  sample  bias  error  is  3imply 


B 


N 


I (A*  - V 

i“l 


(5-3) 


where  N is  the  sample  size.  The  inherent,  statistical  uncertainty  in  the 
sample  bias  is  defined  by  the  confidence  belt  calculated  for  each  sample 
bias  estimate.  If  instead  of  a single  sample  of  size  N,  a large  number  of 
samples  of  size  N are  generated,  then  one  would  expect  the  true  population 
bias  to  fall  within  the  sample  confidence  intervals  for  (1  - a)  percent  of 
the  samples.  The  confidence  intervals  presented  in  these  results  are  com- 
puted from  [18] 


foi  fof~ 

~ - a "W — L_  < B < B + t V — — 

N"1-  N ~ ~ s N-l,_ 2— ’ N 


(5-4) 


where. 


Bg  is  the  sample  bias  computed  in  equation  (5-3), 

B is  the  true,  or  population  bias, 

t i 2_  is  the  students  t-statistic  for  N-l  degrees  of  freedom 
x , *2 

and  a significance  level  of  a 

N is  the  number  of  replications  per  sample  (the  sample  size, 
e.g.  20) 

a is  the  level  of  significance  and  was  chosen  to  be  0.10  for 
this  investigation, 

a?  is  the  sample  variance  of  the  estimated  bearing  as  described 

0 

below. 


■x  ‘’The  simulation  estimates  6^  and  the  theoretical  mean  bearing  esti- 
mate, 0,  described  previously  are  computed  using  the  same  number  of  beam 
positions  for  each  LSE  line  fit. 
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The  sample  variance  of  the  estimated  bearing  is  [8], 


ai 


0 


(5-5) 


I lie  statistical  nature  of  the  sample  variance  or  the  sample  standard  devia- 

i ion,  is  also  reflected  by  its  confidence  interval.  This  confidence 

0 

interval  is  calculated  from  the  following  expression  [8]: 


(K-l) 
1 

N- 1 , 1- 


a 

2 


1 

e 


< 


(5-6) 


vliere  X2  is  the  chi-squared  statistic  for  N-l  degrees  of  freedom 

N-l,  | 

and  a significance  level  of  a, 

N is  the  number  of  estimates  per  sample  and, 
rt  is  the  significance  level,  equal  to  0.10. 

I he  sample  variance  computed  in  equation  (5-5)  reflects  the  random  varia- 
i Ion  of  bearing  estimates  formed  from  p values  of  the  simulation  output  at 
1 corresponding  values  of  the  array  steering  angle  about  the  target  bearing. 

n other  words,  the  variance  in  (5-5)  is  based  on  p observations  of  T 
seconds  of  data.  However,  the  approximate,  theoretical  SBT  bearing 
variance  described  earlier  and  the  Cramer-Rao  variance  bound  derived  in 
chapter  2 are  based  on  only  T seconds  of  observed  data.  To  facilitate 
i direct  comparison  between  the  three  variance  (or  standard  deviation) 
measures,  the  quantity 
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or  * p x o; 
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(5-7) 


is  presented  as  an  approximation  to  the  T-second  bearing  variance  of  the 
ImuJated  SBT  processor.  This  type  of  approximation,  suggested  by 
Overling  [28]  and  adopted  by  Witt  [36],  implies  that  the  angular  region 
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over  which  the  LSE  line  fit  is  made  is  small  enough  to  assume  the  points 
in  the  line  fit  are  p independent  observations  of  statistically  identical 
data. 


5.2  - base! ine  Parameter  Values 

The  results  presented  in  this  chapter  pertain  to  a baseline  set 
of  array,  processor,  signal  and  noise  parameters.  The  baseline  set  of 
parameters  is  listed  in  table  5.1.  The  baseline  values  in  table  5.1  allow 
a realistic  yet  manageable  numerical  analysis  of  the  effects  of  highly 
anisotropic  noise  on  passive  bearing  estimation.  It  is  seen  that  the  only 
variable  parameters  in  the  baseline  study  are  the  number  of  interferers, 
t lie  angular  separation  between  interferers,  and  the  interferer  locations 
relative  to  the  array  axis. 


5.3-  Standard  Deviation  Results 

The  numerical  effects  of  a highly  anisotropic  noise  field  are 
presented  in  figures  5.1  through  5.5  in  terms  of  the  standard  deviation 
of  the  hearing  error  for  the  baseline  case.  Later,  the  bias  error  of 
the  bearing  estimates  is  presented.  The  standard  deviation  of  bearing 
error  or  bearing  estimate  standard  deviation  is  presented  as  a function 
of  the  target-to-interference  separation  which  is  expressed  as  the  phase 
d if  fert-nce. 


where 


Y = (lJJKK  ~ ) * (cos  0T  - cos  Oj)  (5-8) 

^ • (KK-1)  is  the  highest  radian  frequency  in  the 

signal  and  noise  bandwidths  and, 

0j.  is  the  mean  interference  bearing  of  the  cluster  of  1,  2,  or 


A interferers. 


TABLE  5.1 

BASELINE  VALUES  OF  ARRAY,  PROCESSOR, 
SIGNAL,  AND  NOISE  PARAMETERS 


PARAMETER 

SYMBOL 

VALUE 

Number  of  hydrophones 

LL 

10 

Hydrophone  separation 

D 

20  feet 

Integration  time 

T 

0.128  seconds 

Sample  rate 

250  Hz 

Sign.i  l-to--Noise  ratio 

SNR 

1.0 

True  target  bearing 

°T 

90° 

Inter ference-to-Noise 
ratio 

1NR 

0.1 

Number  of  interferers 

MM 

1,  2,  4 

Interferer  separation 

(A6) 

3,  6,  9 degrees 

Interference  location 

9I 

Theoretical:  0-180  degrees 
Simulation:  30,  45,  60,  70 
76,  82,  90  degrees 

Signal  and  noise  process 
bandwidth 

101.5  Hz 

Discrete  bandwidth 

KK-1 

13 

Number  of  beam  positions 
per  estimate 

ISTEER 

Cp) 

13 

Angular  steering  sector 
covered 

86  - 94  degrees 
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Since  l lie  standard  deviation  of  error  metrics  are  symmetric  about  the 
zero  phase  separation,  these  results  are  presented  as  function  of  |y| 
rattier  than  Y. 

Figure  5.1  shows  the  standard  deviation  of  error  for  the  case  of 
a single  plane  wave  interference.  Figures  5.2a,  b,  and  c show  the  results 
for  two  interferers  separated  by  3,  6,  and  9 degrees,  respectively.  Like- 
wise, figures  5.3a,  b,  and  c show  the  results  for  four  interferers.  Three 
general  observations  are  readily  apparent  from  these  results.  First,  in 
every  case  the  standard  deviation  of  error  exhibits  a different  behavior 
In  each  of  three  segments  of  the  range  of  |y|  values.  For  remote  inter- 
feres, all  three  error  measures  are  independent  of  the  target- to-inter- 
ference  separation  | Y | , and  only  somewhat  larger  than  the  standard  devia- 
tion for  isotropic  noise  alone.  For  intermediate  values  of  |y|  (approxi- 
mately 0.4  to  1.0),  the  standard  deviation  metrics  exhibit  a peak  value 
(or  a relative  maxima).  Finally  at  small  separations  of  target  and  inter- 
fe.rer,  the  behavior  of  the  three  metrics  appears  to  vary  from  case  to  case 
as  a complex  function  of  MM,  A0  and  |y| . The  detailed  behavior  of  the 
standard  deviation  of  error  in  these  three  segments  will  be  examined  later. 

A second  general  observation  is  that  the  Cramer-Rao  bound  indeed 
provides  a lower  bound  on  the  bearing  error  except  at  very  small  phase 
separations  (small  |y|),  thus  illustrating  the  performance  degradation  of 
the  SBT  processor  over  the  ML  processor.  In  some  cases,  at  small  |y|  values 
the  theoretical  SBT  error  becomes  less  than  the  Cramer-Rao  bound.  In 
these  regions,  the  target  and  one  or  more  interferers  are  unresolvable  by 
both  Ml.  and  SBT  processors.  However,  the  SBT  processor  essentially  "sees" 
a single  source  of  directional  noise  (with  increased  level)  and  a background 
noise  field  that  is  nearly  incoherent,  or  isotropic.  As  a result  of  the 
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Figure  5.1  - Standard  Deviation  of  Estimation  Error,  MM=1 
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Figure  5.2  - Standard  Deviation  of  Estimation 
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higher  apparent  signal  level,  it  is  seen  in  figure  5.1  that  the  SBT 
error  al  small  |y|  falls  below  its  own  isotropic-noise  error  level.  The 
performance  of  the  ML  estimator  which  is  bounded  by  the  larger,  CR  bound 
reflects  the  tacit  assumption  that  the  optimal  estimator  is  perfectly 
matched  to,  or  has  perfect  knowledge  of,  the  true  background  noise  field. 
Hence,  the  fact  that  the  SBT  and  ML  processors  are  addressing  "different" 
noise  fields  at  small  |Y|  values  is  reflected  by  the  seemingly  superior 
performance  of  the  SBT  processor. 

The  third  general  observation  is  that  the  bearing  error  of  the 
simulated  SBT  is,  at  the  chosen  confidence  level  of  0.90,  significantly 
poorer  than  the  approximate,  theoretical  SBT  bearing  error.  This  appears 
true  in  each  of  the  figures  5.1  through  5.3  at  most  target-to- interference 
separations.  The  appearance  that  this  difference  is  almost  constant  over 
all  | Y | values  is  explained  by  the  use  of  the  same  set  of  input  normal 
deviates  for  each  of  the  seven  target-to- interference  separations  examined 
by  the  simulation.  The  difference  between  the  simulation  and  theoretical 
SBT  bearing  errors  occurs  consistently  throughout  the  investigation. 
Although  the  exact  cause  of  this  somewhat  unexpected  difference  is  unknown 
a combination  of  one  or  more  of  the  following  points  is  suggested  as  a 
possible  explanation: 

1)  The  number  of  replications  (i.e.,  simulation  runs)  is  small  rela- 

tive to  the  number  of  points  used  in  the  line  fit.  The  SBT  output 
values  used  in  the  LSE  line  fit  occur  in  an  eight  degree  segment  of 
steering  angles  between  86®  and  94®.  The  true  bearing  error,  which  is 
estimated  by  the  simulation  results,  is  on  the  order  of  one  to  two 
degrees  - a significant  fraction  of  the  8®  segment  of  steering  angles 
forming  the  line  fit.  Hence,  the  number  of  points  used  to  form  each 


estimate  and  the  number  of  estimates  per  sample  will  signif icar tiy 
affect  the  standard  deviation  of  estimation  error. 

2)  The  theoretical  assumption  that  the  SBT  output  is  a linear  funcr ;c 
of  steering  angle  is  Invalid  for  the  86-94  degree  range  over  which  the 
simulation  forms  the  LSE  straight  line  fit.  Figure  5.4  shows  a plot  S 

the  theoretical  SBT  output  mean,  z , and  standard  deviation,  a , as 

o 

functions  of  steering  angle,  a,  using  the  expressions  derived  in 
chapter  3.  Notice  that  in  the  86-88  and  92-94  degree  regions,  the 
assumption  of  linearity  is  somewhat  strained.  The  nature  of  this  non- 
linearity would  tend  to  cause  a lower  LSE  line  slope  than  expected  and 
thus  a higher  standard  deviation  of  error. 

3)  There  are  13  points  used  by  the  LSE  line  fit  estimator  to  obtain 
the  baseline  simulation  bearing  estimates.  Under  the  assumption  that 
13  T-second  observations  of  statistically  identical  random  data  produce 
an  error  variance  which  is  13  times  lower  than  that  obtained  for  a 
single  T-second  observation,  the  standard  deviation  of  error  has  been 
multiplied  by  /l3  for  presentation  in  figures  5,1  to  5.3.  Here  lies 
the  third  possible  contributing  cause  of  the  difference  between  simu- 
lated and  theoretical  bearing  error.  That  is,  the  observations  at 

13  different  beam  positions  (86°  - 94°)  may  not  be  statistically 
identical,  at  least  to  the  extent  that  it  is  optimistic  to  assume  that 
the  simulation  error  variance  should  be  13  times  lower  than  an  "equiva- 
lent" T-second  error  variance.  A less  optimistic  assumption  would,  of 
course,  result  in  lower  "T-second  equivalent"  simulation  bearing  errors 
than  those  presented  in  figures  5.1-5. 3. 
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4.'  The  actual  observation  time  of  0.128  seconds  is  too  small  for  the 

approximation,  made  in  the  theoretical  SBT  standard  deviation  o?  erro- 
derivation,  that  the  spectral  density,  X(w),  can  be  expressed  as 

X (u)  * ~~  ' E jx(jw)  x (jw)  j (5-9) 

when-  x(ju»)  is  the  Fourier  Transform  of  the  signal-plus-noise-plus- 
inteif erence  process  existing  in  the  water  at  the  face  of  the  array. 

Recalling  the  first  general  observation  drawn  from  the  standard 
deviation  of  error  results,  it  seems  appropriate  to  discuss  the  detailed 
eifects  of  a highly  anisotropic  noise  field  in  terms  of  the  three  levels  of 
target- l.o- interference  separations  — remote  interferers,  intermediate  or 
"peak"  interferers,  and  nearby  or  coincident  interferers.  In  the  inter- 
mediate | Y J - region,  the  location  of  the  peak  error  varies  somewhat  with 
and  fO.  Also,  the  SBT  error  and  the  CR  bounds  peak  at  different  loca- 
tions. For  the  single  interferer  shown  in  figure  5.1,  the  locations  of  the 
ei ror  peaks  are  identified  with  the  resolution  beamwidths  of  the  SBT  and 
Mi  processors.  In  a previous  theoretical  investigation  [7],  Cox  derived 
and  evaluated  an  expression  for  the  resolution  power  of  the  optimal 
(maximum  likelihood)  beamformer.  From  his  results  ([7],  figure  4,  j ...  781) 
it  is  straightforward  to  determine  that  the  resolution  "beamwidths"  (i.e. 
ri solution  power  in  terms  of  angular  separation)  for  the  parameters  of 
interest  are  14.2°  for  a conventional  beamformer  and  10°  for  the  optimal 
beamformer.  From  figure  5.1,  we  see  that  the  peak  bearing  error  occurs 
when  the  target  and  interferer  are  separated  by  a resolution  beanjwidth. 

The  results  in  figures  2 and  3 indicate  that  for  multiple,  closely  spaced 
(e.g.  3°)  interferers,  the  standard  deviation  of  error  peaks  when  the 
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mean  Interference  location  is  approximately  one  resolution  beamwidth  from 
the  target.  As  the  interference  separation,  A0,  increases  the  | Y j posi- 
tion of  the  peak  error  increases,  reflecting  the  transition  to  indepena  nr 

interferer  effects  at  A0  values  greater  than  a resolution  beamwidth. 

The  effects  of  increasing  the  number  and  mutual  separation  of 

remotely  positioned  interfering  noise  sources  is  summarized  in  fignra  5.5. 

It  Is  obvious  here  that  the  mutual  separation  of  interferers,  like  their 
position  or  bearing,  has  no  effect  on  the  standard  deviation  of  error 
when  the  target- to- interference  separation  is  large.  The  CR  bound  indi- 
cates that  the  ML  estimator  will  also  be  independent  of  the  number  of 
remote  interferers.  The  simulation  and  theoretical  SBT  errors  indicate, 
however,  that  the  sub-optimal  processor  performance  worsens  with  increases 


in  the  number  of  remote  interferers.  Since  the  SBT  is  designed  to  discrim 
inate  against  incoherent  noise,  it  cannot  cancel  or  "de-emphasize"  the 
remote  interference  which  in  turn  increases  the  equivalent  background 
noise  level.  As  discussed  previously,  the  simulated  SBT  bearing  error  is 
significantly  higher  than  the  theoretical  SBT  error.  However,  its  depend- 
ence on  MM  and  A0  closely  mirrors  that  of  the  theoretical  measure. 

The  behavior  of  the  standard  deviation  of  error  in  the  mid  or 
intermediate  range  of  |y|  values  is  shown  in  figure  5.6.  This  expresses 
the  peak  error  as  a function  of  A0.  As  one  might  expect,  the  peak  values 
of  both  the  CR  bound  and  the  SBT  error  for  multiple  interferers  approaches 
the  single-interference  values  as  the  mutual  separation  increases.  For 
interference  separations  on  the  order  of  two  times  the  resolution  beam- 
width,  multiple  and  single  interferers  would  yield  the  same  peak  values. 
The  CR  bound  now  shows  a slight  dependence  on  number  interferers  while  the 
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SBT  processor  shows  a high  dependence  on  number  of  interferers  when  tbi:. 
separat  ion  is  small. 

The  difference  in  the  slop's  of  the  error  curves  in  figure  5.6 
Indicates  a rapid  deterioration  in  the  relative  performance  (i.e.  the  rat.iv. 

1 CUT 

-)  of  the  SBT  processor  as  the  number  of  closely  spaced  interferers 

°CR 

increases.  MacDonald  and  Schultheiss  [19]  have  shown  that  for  incoherent 
background  noise  (no  interferers)  the  relative,  theoretical  degradation 
of  the  SBT  processor  is 

CR  v ' 

where  l.L  is  the  number  of  uniformly  spaced  hydrophones.  For  the  10- 
hydrophone  array  of  interest,  this  ratio  is  1.15.  The  results  in  figure 
5.1  show  theoretical  and  simulation  degradations  relative  to  the  CR  bound 
for  a single  interferer  of  1.23  and  1.47,  respectively.  For  four  inter- 
ferers separated  by  three  degrees,  the  relative  performance  degradation 
has  increased  to  1.69  and  2.39  for  the  theoretical  and  simulation  measures 
respectively.  For  smaller  separations,  we  would  expect  larger  relative 
performance  losses.  This  data  represents  a very  serious  degradation  in 
bearing  estimation  performance  of  the  sub-optimal  processor  over  the 
optimal  processor  but  is  not  unexpected  in  light  of  the  sub-optimal  pro- 
cessor design  emphasis  on  isotropic  noise  fields. 

Figure  5.7  shows  the  standard  deviation  of  bearing  error  behavior 
for  coincident  target  and  interference.  As  noted  earlier,  the  CR  bound 
appears  not  to  provide  a lower  bound  on  the  bearing  estimation  error  for 
all  separations  of  multiple  interferers.  Since  the  SBT  processor  has  no 
knowledge  of  any  directional  aspects  of  the  noise  field,  an  interference 
in  the  proximity  of  the  target  is  interpreted  as  an  increase  in  signal 
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Figure  5.7  - Standard  Deviation  of  Error  - Coincident  Interference 
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Level.  Consequently , the  SRT  processor  sees  a target  with  an  effective 
signal  level  higher  than  the  true  target  level. 

Since  Jie  target  is  coincident  with  the  mean  bearing  of  the  inter- 
ference duster,  the  standard  deviation  of  error  is  seen  to  increase  with 
Increasing  A0  for  multiple  interferers.  One  would  expect  this  beha/ior 
for  A6  less  than  or  equal  to  the  resolution  beamwidths.  As  A0  increases 
above  a resolution  beamwidth,  there  can  be  at  most  one  interference  within 
the  beam  and  the  effects  of  the  remaining  members  of  the  interference 
cluster  will  decrease.  Also,  the  simulation  standard  deviation  of  error 
is  again  seen  to  mirror  the  theoretical  metric  but  at  a somewhat  higher 
level. 

5.4  - Bias  Error  Results 

The  mean,  or  bias,  errors  determined  for  the  baseline  study  are 
presented  in  figure  5.8  for  a single  interferer  and  figures  5.9  and  5.10 
(a,  b,  and  c)  for  two  and  four  interferers  at  3,  6,  and  9 degree  separa- 
tions, respectively.  The  bias  curves  in  figures  5.8  - 5.10  reflect  only 
SBT  errors  since  the  optimal  (ML)  estimator  is  by  assumption  unbiased. 

Several  general  observations  can  be  made  by  examination  of  figures 
5.8  - 5.10. 

1)  The  bias  error  is  anti-symmetric  about  0^,  and  of  significant 
magnitude  over  a wide  range  of  target-to- interference  separations  00 . 
As  one  might  expect,  the  estimated  bearings  are  biased  in  the  direction 
of  the  interference  location. 

2)  The  bias  errors  exhibit  peaks  (relative  extrema),  whose  locations 
and  widths  vary  noticeably  with  number  and  separation  of  interferers. 
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For  multiple  interferers,  the  peak  values  decrease  with  increasing 
interference  separation  while  the  peak  widths  increase  with  increasing 
separation. 

3)  The  simulation  and  theoretical  bias  errors  are  in  close  agreement 

for  all  cases.  Unfortunately,  this  agreement  offers  no  solid  clue  as 
to  the  source  of  discrepancy  between  the  theoretical  and  simulated 
SBT  standard  deviations  of  error. 

The  detailed  behavior  of  the  absolute  value  of  the  peak  SBT  bias 
error  as  a function  of  number  and  separation  of  interferers  is  shown  in 
figure  5.11.  The  close  agreement  between  the  theoretical  and  simulation 
bias  metrics  is  apparent  in  figure  5.11.  The  bias  measures  decrease 
linearly  as  A 9 increases.  The  rapid  performance  degradation  of  the  SBT 
processor  as  the  number  of  closely  spaced  interferers  increases  is  evi- 
denced by  the  larger  negative  slope  of  the  MM  * 4 bias  curve.  A similar 
behavior  was  witnessed  for  the  peak  standard  deviation  of  error  curves 
presented  in  figure  5.6.  As  with  the  peak  standard  deviation  of  error, 
the  peak  bias  error  for  multiple  interferers  will  approach  that  for  a 
single  interferer  as  A0  becomes  large. 


CHAPTER  Vi 


NUMERICAL  RESULTS  - PARAMETER  EXTENSIONS 

Tin-  results  presented  in  chapter  5 provide  basic  insight  into  the 
efforts  of  a highly  anisotropic  noise  field  on  bearing  estimation  by  optimal 
(ML)  and  sub-optimal  (SBT)  processors.  In  this  chapter,  added  insight  into 
the  effects  of  such  a noise  field  is  provided  by  the  numerical  results  of 
ox  Lending  the  range  of  values  of  selected,  fixed  baseline  parameters.  To 
maintain  manageability  and  economic  feasibility,  only  two  of  the  potentially 
more  interesting  parameters  are  examined  in  this  phase  of  the  study.  The 
dependency  of  hearing  estimation  performance  on  received  signal-to-noise 
(SNR)  ratio  is  presented  in  section  6.1  with  the  numerical  results  for  SNR 
values  of  0.5  and  2.0.  1’he  following  section  presents  the  results  of  ex- 
tending the  array  size  to  20  and  40  hydrophones  (for  an  SNR  of  1.0).  Except 
as  noted  in  each  section,  the  remaining  fixed  and  variable  parameters  assume 
their  baseline  values  as  defined  in  table  5.1. 

(>.  1 - Extension  of  Signal-to-Noise  Ratio 

Ibis  section  presents  the  results  of  extending  signal-to-noise 
tange  i<>  0.5  and  2.0.  This  SNR  range  covers  many  of  the  situations  en- 
countered in  aetuaL  target  tracking  applications.  Although  the  simulated 
hearing  error  has  been  consistently  higher  than  the  approximate  theoretical 
error,  its  overall  behavior  as  a function  of  number  and  location  of  inter- 
ferons is  again  expected  to  mirror  the  theoretical  metric.  Consequently, 
tin'  simulation  was  not  exercised  for  Interference  bearings  of  30°  and  60°. 
Tlio  simulation  results  for  the  remaining  five  interferer  positions  (45°, 
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70",  /(>”,  82°,  90°)  are  presented  in  this  section.  Furthermore,  at  each 
interlerer  position  the  simulation  sample  size  is  limited  to  20  replica- 
lions.  N<  lliuL  throughout  the  SNK  parameter  extension,  the  array  size 
remains  at  10  hydrophones. 

The  numerical  results  o£  evaluating  the  three  standard  deviation 
ol  error  meLrics  and  the  two  bias  error  metrics  are  presented  in  figures 
1-7  and  8-14,  respectively  of  appendix  D.  For  clarity,  the  baseline  resul 
al  SNK  = 1.0  are  not  repeated  in  these  figures. 

Several  observations  can  be  made  upon  examination  of  the  standard 
deviation  of  error  metrics  in  figures  D.l  - 0.7  regarding  the  extended 
SNK  parameter  values: 

1)  Over  the  entire  range  of  |y|,  the  standard  deviation  of  error  in- 
creases  at  small  s Lgnal-to-noise  ratios.  For  the  most  part,  this 
phenomena  results  from  the  loss  in  signal  power  relative  to  the  isotrop 
i<-  component  of  the  background  noise  level.  Another  contribution  to 
this  effect  is  possible  since,  even  at  remote  interference  bearings, 
the  total  Interference  power  will  contribute  somewhat  to  the  received 
background  noise  level  seen  by  the  two  processors. 

2)  The  larger  peak  magnitude  of  the  standard  deviation  of  error 

and  the  seemingly  superior  coincident  (|y|  = 0)  performance,  bv  the  SBT 
processor  relative  to  the  OR  bound  at  low  SNK,  results  from  a loss  in 
signal  level  relative  to  the  received  interference  power.  In  other 
word:;,  these  are  phenomena  observed  as  a result  of  a decreased  nignal- 
to  i ol  ••r|  ereiice  power  ratio.  Observe  also  that,  although  the  error 
peaks  are  sharper  (larger  relative  t.o  the  bound)  at  SNR  = 0.5,  their 
he.ning  location  ( | Y | value)  remains  constant  over  the  SNR  range  of 
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3)  The  apparent  discrepancy  between  the  theoretical  and  the  simulated 

SBT  standard  deviations  of  error  chat  was  discussed  in  chapter  5 is 
also  maiilest  in  the  results  presented  in  appendix  D.  It  is  noteworthy, 
however,  that  the  discrepancy  worsens  significantly  around  the  error 
peaks  for  multiple  interferers,  particularly  when  closely  spaced.  For 
example,  consider  figure  D.5  for  four  interferers  separated  by  3°. 

Tlu-  simulation  errors  at  |y|  equal  0.355  and  0.618  for  SNR  - 0.5  h3ve 
been  adjusted  to  reflect  smaller  sample  sizes  of  18  and  16  replications, 
respectively.  These  sizes  result  from  an  after-the-fact  rejection  of 
three  obviously  unacceptable  bearing  estimates.  The  small  "x"  at  Y 
equal  0.355  and  0.618  indicate  that  respective  error  magnitudes  of 
4.8  and  12.24  degrees  were  obtained  for  the  original  simulation  sample 
sizes.  The  rationale  here  is  that  such  poor  bearing  estimates  (6  = 81.9°, 
78. 6°,  and  75.2°)  would  be  ignored  in  any  automatically  or  manually 
operated  bearing  estimation  system.  Furthermore,  it  is  noted  that  the 
full,  original  sample  sizes  of  19  and  18  replications  at  |y|  equal  to 
0.355  and  0.618,  respectively,  and  SNR  = 0.5  in  figure  D.5  are  smaller 
than  the  other  20-run  samples  due  to  an  absence  of  a zero-crossing  by 
the  simulation  output,  and  hence  no  bearing  estimate,  for  three  simula- 
tion runs. 

The  unacceptable  estimates  from  some  runs  and  the  absence  of  a zero- 
crossing  in  other  runs  is  a consequence  of  the  simulated  SBT  output  be- 
coming non-linear  in  the  86°-94°  region  of  the  T.SE  line  fit.  This  non- 


linearity is  also  responsible  for  the  larger  discrepancy  between  the 
theoretical  and  simulated  standard  deviations  of  error  at  these  |y|  - 
values.  The  nature  of  this  non-linearity  is  shown  in  figure  6.1.  This 
figure  shows  the  theoretical  mean  and  standard  deviation  of  the  SBT 
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output  (equal  ions  (3-22)  and  (3-23))  for  an  SNR  of  0.3  and  four  inter- 
fotors  separated  by  3°  and  centered  at  a bearing  of  76°  (|y|  = 0.618). 

Tli.  occurrence  of  the  non-linearity  in  figure  6.1  for  small  zq  values 
indicates  not  only  a large  bearing  estimation  variance  but  also  an 
asymmetric  probability  density  function  (pdf)  for  0^.  The  unacceptably 
lov,.  bearing  estimates  mentioned  previously  (81.9°,  78.6°,  and  75.2°)  are 
a consequence  of  the  unusually  long  "tail"  in  the  0^  pdf.  As  seen  in 
figure  6.1,  the  lower  one-sigma  value  of  the  0^  pdf  occurs  approximately 
8./"  below  the  86.7°  theoretical  mean  while  the  upper  one-sigma  value 
occurs  only  3.5°  above  the  mean.  Hence,  while  the  three  exceptionally 
low  hearing  estimates  might  be  unacceptable  in  an  operational  applica- 
tion, they  are  statistically  valid  estimates  which  are  indicative  of 
a lower  SNR  limit  on  bearing  estimation  by  the  SBT  processor. 

Several  general  observations  can  also  be  made  from  the  SBT  bias 

ci roi  results  presented  in  figures  D.8  through  D.14: 

I?  SNR  variation  affects  the  SBT  bias  error  primarily  by  compressing 

or  .ciipl  lfying  the  bias  error  peaks.  At  low  SNR  levels  (approximately 
equal  to  the  interference  level)  the  bias  peak  has  risen  to  the  extent 
indicative  of  bearing  estimates  midway  between  target  and  interference. 
Thi  is  approximately  true  for  SNR  « 0.5  in  figure  D.12.  For  an  SNR 
lev  than  the  INR  at  close  separations,  the  SBT  processor  would  be,  in 
os; .ice,  estimating  the  interference  bearing  rather  than  the  signal 
lie.  , i i ng. 

2)  For  low  bias  values  such  as  at  SNR  =•  2.0  and  MM  1,  a difference  in 

tlu  theoretical  and  simulation  bias  metrics  appears.  However,  at  a 90% 
coni  ulence  level  we  cannot  say  that  this  difference  is  significant.  Also 
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noli:  Lli.it  lor  MM  - 4 anil  SNR  - 2.0  (figure  0.12)  the  theoretical  and 
simulation  bias  metrics  tend  to  converge  as  their  magnitudes  increase. 
'Mies.*  i w<  cents  mighL  suggest  the  difference  in  bias  metrics  at  high 
SNR  is  statistical  (i.e.  due  to  a small  sample  size  and/or  too  few 
points  in  LSK  line  fit). 

i)  Recalling  the  non-linearity  of  the  mean,  zq  and  standard  deviation, 

0 , of  the  SBT  output  shown  in  figure  6.1,  one  would  expect  it  to 
o 

cvenl.ua l. Ly  affect  the  bearing  estimation  bias  error  as  it  did  the  stand- 
ard deviation  of  bearing  error.  Indeed,  one  observes  in  figure  D.12 
ih.it  the  absolute  values  of  simulation  bias  error  are  greater  than  the 
theoretical  values  for  the  uncorrected  sample  sizes  (19  and  18  replica- 
tions) at  V = i 0.355  and  ± 0.618.  Removal  of  the  three  unacceptable 
ci  iinaLos  as  discussed  previously  yields  simulation  bias  errors  somewhat 
closer  to  the  theoretical  metric.  From  figures  D.5  through  D.12,  it  is 
seen  that  for  four  closely  spaced  interferers  located  within  about  20° 

01  a 0.5  SNR  target  the  SBT  processor  is  practically  useless  as  a 
bearing  estimator. 

The  behavior  of  the  standard  deviation  metrics  as  functions  of 
sign  i l-to-noise  ratio  can  also  be  examined  parametrically  in  each  of  the 
i lire  previously  defined  regions  of  target-to- interference  separation,  ]y|. 
figure  6.2  shows  the  effects  of  increasing  the  signal-to-noise  ratio  for 
remote  interference  bearings  (large  |y|).  The  primary  observation  from 
figine  6.2  is  that,  except  for  the  simulation  error  at  SNR  * 2.0,  the 
standard  deviation  of  error  metrics  vary  exponentially  with  SNR.  For 
lull  logarithmic  scaled  coordinates,  it  is  easy  to  show  that  the  unsealed 
slop  of  a straight  line  is  the  exponential  power  of  the  abscissa  which 
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governs  the  variation  of  the  ordinate  function.  Consequently,  we  see  from 
I iguto  6.2  the  standard  deviation  of  error  varies  inversely  with  the  square 
loot  o!  the  fractional  signal-to-noise  ratio,  SNR.  This  contrasts  with 
MacDonald's  [18|  results  that  the  standard  deviation  of  error  (remote 
interference  bearings)  is  inversely  proportional  to  SNR  for  conditions  of 
high  INK  and  low  SNR. 

The  CR  bound  is  seen  to  be  identical  for  all  values  of  MM  and 
AO  at  remote  interference  bearings.  The  SBT  metrics,  while  also  independ- 
ent of  AO,  show  a slight  dependence  on  the  number  of  interferers  present,  MM. 

The  simulation  error  metric  in  figure  6.2  appears  approximately 
constant  over  the  1.0  to  2.0  range  of  SNR  values,  thus  indicating  a larger 
l heoret leal /simulation  error  discrepancy  than  previously  identified.  This 
unexpected  behavior  is  probably  indicative  of  a bearing  resolution  limit 
characteristic  of  our  digital  estimation  algorithm.  In  other  words,  the 
simulation  random  errors  at  SNR  = 2.0  based  on  13  T-second  observations 
(13  beam  positions)  are  not  resolvable  beyond  one-half  the  angular  separa- 
I ion  between  adjacent  beam  positions  (0.33  degrees).  Consequently,  the 
estimation  algorithm  would  yield  erroneous  results  when  the  true  random 
error  falls  below  half  the  beam  position  spacing.  The  relatively  small 
sample  size  at  SNR  = 2.0  (20  replications)  would  also  contribute  to  this 
behavior.  If  these  theoretical/siraulation  discrepancies  are  indeed  caused 
by  the  LSE  line  fit,  then  hopefully  they  can  be  eliminated  by  proper 
select  ion  of  the  parameters  of  the  bearing  estimator.  It  is  clear  that 
I here  is  room  for  more  investigation  into  the  digital  implementation  of 
bearing  estimators  which  are  based  on  zero  crossings  of  an  array  processor 
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output . 
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Figaro  6.3  indicates  that  the  peak  standard  deviation  of  error  is 
1 1 so  inversely  proportional  to  the  square  root  of  the  signal-to-noise  ratio, 
unly  i no  3°  separation  of  multiple  interferers  is  shown  in  figure  6.3.  The 
peak  errors  for  6 and  9 degree  separations  are  smaller  than  the  40  = 3° 
result  i hut  exhibit  similar  behavior  as  SNR  increases.  The  SBT  metrics  in 
ilgui<  6.3  indicate  a greater  sub-optimal  processor  degradation  as  the 
numb  i of  closely  spaced  interferers  increases  than  indicated  by  the  CR 
hound  tor  the  Ml.  estimator.  The  simulated  SBT  error  appears  to  behave  more 
like  i he  theoretical  SBT  error  at  the  higher,  peak  error  values  than  it 
did  it>i  remote  interference  bearings. 

For  the  case  of  a coincident  target  and  interference,  the  behavior 
t>f  the  standard  deviation  of  error  metrics  shown  in  figure  6.4  closely 
tosemhlos  that  of  the  error  metrics  for  remotely  positioned  interferers 
(figure  6.2).  Notable  differences  seen  in  figure  6.4  are: 

1)  Generally  lower  (by  0.3°  to  0.5°)  SBT  standard  deviation  of 

error  metrics. 

2)  The  CR  bound  depends  on  the  number  of  interferers. 

3)  The  SBT  error  is  not  consistently  bounded  below  by  the  CR 

bound . 

As  a final  measure  of  the  SNR  dependence  of  bearing  estimation  per- 
lorm.une  metrics,  consider  the  peak  SBT  bias  errors  presented  in  figure 
(>.5.  The  previously  noted  theoretical/simulation  bias  discrepancy  is  again 
> v i di  n i at  SNR  = 2.0.  This  bias  discrepancy  decreases  as  the  magnitudes  of 
«ach  m.  trie  increase.  This  discrepancy  was  previously  associated  with  a 
hoar  in)1,  resolution  limit  of  the  SBT  simulation  and  the  small  sample  size 
it  SNi:  = 2.0.  Of  primary  importance,  however,  are  the  following  two 
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Figure  6.3  - Peak  Standard  Deviation  of  Error  vs.  SNR 
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1)  Kor  I,  2,  and  4 closely  spaced  Interferers,  the  SBT  bias  error 

is  inversely  proportional  to  Lhe  fractional  signal-to- 
noise  ratio,  SNR. 

2)  The  SBT  bias  errors  are  directly  proportional  to  the  number,  MM, 

of  closely  spaced  interferers. 


These  iwo  observations  indicate  that  the  peak  SBT  bias  error  is,  to  a good 

approximation,  inversely  proportional  to  the  signal -to-total-coherent- 

no,se  level,  • 

N 

6.  ' - Kxtended  Array  Size 

This  section  presents  the  results  of  extending  the  array  size,  LL, 
to  20  and  40  hydrophones.  These  results  address  only  the  case  of  a signal- 


lo  noise  ratio  of  1.0.  Several  adjustments  to  the  baseline  set  of  parameter 


values  wore  necessitated  by  the  extension  of  the  array  size.  These  adjust- 
ments are  required  by  the  smaller  resolution  beamwidths  obtained  with  the 
larger  apertures.  For  example,  since  adjacent,  resolvable  Interferers  will 
all  eel  i lie  array  processors  almost  independently,  it  is  of  little  interest 
to  study  multiple,  widely  separated  (much  more  than  a bearawidth)  interferers 
and  since  Lite  larger  20-  and  40-hydrophone  arrays  have  beamwidths  smaller 
thin  <1°,  the  interference  separation  parameter  assumes  only  the  3 and  6 
degree  values  for  this  portion  of  the  analysis. 

A related  aspect  of  the  larger  array  size  is  the  smaller  segment 
of  steering  angles  over  which  the  SBT  output  is  approximately  linear. 

Figures  b.6  and  6.7  show  the  approximate  theoretical  output  mean  and  stand- 
ard deviation  for  array  sizes  of  20  and  40  hydrophones,  respectively.  In- 
creasing the  array  size  appears  to  amplify  the  SBT  output  extrema  and  shift 
them  closer  to  the  zero-crossing  and  hence  decrease  the  angular  region  of 

II 
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linearity.  Kor  this  reason,  the  simulation  estimates  were  formed  from 
output  values  in  the  88-92  degree  region  for  LL  = 20  and  in  the  89-91  degree 
region  for  LL  = 40. 

The  extreme  costs  associated  with  exercising  the  SBT  simulation  for 
large  arrays  mandate  still  other  parameter  changes.  The  cost  of  exercising 
a digital  computer  simulation  results  mainly  from  processor  (CPI1)  run  time 
and  core  memory  usage.  Figure  6.8  illustrates  the  time  and  memory  require- 
ments lor  exercising  the  simulation  on  Catholic  University's  DEC  KL-10 
computer.  Processing  time  in  CPU  minutes  is  given  by  the  left  ordinate 
while  memory  usage  in  kilo-core-seconds  is  given  by  the  right  ordinate. 

A k i lo-eore-second  is  defined  as  1024  words  of  CPU  memory  used  per  second 
of  run  time.  Figure  6.8  reflects  the  costs  required  to  obtain  20  bearing 
estimates,  at  13  beam  positions  per  estimate  for  a single  value  of  the  mean 
Inter  I erenee  bearing,  8 . The  core  memory  usage  and  CPU  time  are  respec- 
tively about  100  and  9.5  times  greater  for  40  hydrophones  than  for  10  hydro- 
phone?;. Due  to  these  higher  costs,  only  nine  beam  positions  are  used  to 
form  each  of  the  20-  and  40-hydrophone  bearing  estimates.  In  this  investi- 
gation, we  must,  unfortunately,  accept  the  expected  increase  In  bearing 
error  associated  with  a decrease  in  the  number  of  points  in  the  LSE  line 
fit.  In  subsequent  Investigations,  it  might  be  possible  to  circumvent  this 
trade  off  by  optimizing  the  simulation  software  programs  and  by  even 
employing  special,  dedicated  hardware  in  some  applications. 

Other  required  parameter  changes  are  the  mutual  interference 
separai ions  and  target-to-interf erenee  separations  that  were  addressed  by 
the  simulation  exercises.  Unlike  the  theoretical  metrics,  the  simulation 
error  metrics  are  only  evaluated  for  mutual  interference  separations  of 
3"  and  also  for  only  two  separations  of  target  and  interference  (i.e.  ]y|). 
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These  wore  the  "peak  error"  Interference  separations  and  the  case  of 
coincident  target  and  Interference.  Recall  that  in  the  baseline  investiga- 
tion results  for  remotely  positioned  interferers  (figure  5.5)  the  SBT 
standard  deviation  of  estimation  error  is  independent  of  A0  and  has  a 
small  linear  dependence  on  MM.  Since  similar  results  are  expected  for  the 
extended  array  sizes,  the  bearing  error  for  remotely  positioned  interferers 
was  only  obtained  for  the  MM  = 1 case. 

To  summarize,  the  theoretical  results  are  presented  for  A0  equal 
'1°  and  6°  while  the  simulation  results  are  presented  only  for  A0  * 3°. 
Furthermore,  each  simulation  estimate  for  20  and  40  hydrophones  is  formed 
from  the  SBT  outputs  at  nine  beam  positions  spanning,  respectively,  the 
88-92  and  89-91  degree  sectors  of  steering  angles.  The  sample  size  remains 
at  20  runs  per  sample. 

The  estimation  standard  deviation  and  bias  metrics  are  presented  in 
figures  1-5  and  figures  6-10  of  appendix  E,  respectively.  In  each  figure, 
the  baseline  results  for  a 10-hydrophone  array  are  repeated  for  comparison 
with  the  20-  and  40-hydrophone  results.  An  immediate  observation  from 
these  figures  is  the  profound  effect  that  the  decrease  in  resolution  beam- 
width  has  on  both  standard  deviation  and  bias  error  metrics.  These  effects 
are  the  shift  of  the  error  peaks  in  the  direction  of  the  target  bearing  and 
the  "(isc 1 1 latory"  behavior  of  the  metrics  for  multiple,  resolvable  inter- 
ferers. From  Cox's  previous  study  of  target  resolution  [7],  the  optimal 
and  conventional  bearawidths  have  been  computed  for  the  parameters  of 
interest  and  are  shown  in  table  6.1.  This  table  and  the  results  in  appendix 
K indicate  that  the  standard  deviation  of  error  peaks  for  single  and  multiple, 
closely  spaced  interferers  occur  approximately  one  resolution  bearawidth 
from  the  target  bearing.  Widely  spaced  interferers  should  affect  the  ML 
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TABLE  6.1 

RESOLUTION  BEAMWIDTHS  FOR  EXTENDED 
ARRAY  LENGTHS  [7] 


Beamformer 

Resolution 

Phase 

Conf iguration 

Beamwidth 

Separation,  j 

y| 

(//  Sensors) 

10 

O 

O 

rH 

0.44 

Opt imal 

20 

4.6° 

0.20 

Beam former 

40 

1.96° 

0.09 

10 

14.2° 

0.63 

Convent ional 
Beamformer 

20 

7.1° 

0.32 

40 

3.55° 

0.16 

’The  tabulated  |y|  is  the  phase  separation  between  the  target 


< il  0.|t  - 90°)  and  an  interference  cluster  centered  at  one  resolution  beam- 


width  in  cither  side  of  it. 


. ml  SBT  processors  independently  and  result  in  multiple,  equally  spaced 
I oaks  in  tlie  standard  deviation  of  error  metrics.  This  "picket  fence" 

< I'fect  is  indeed  observed  in  figure  E.5  for  LL  = 40.  When  each  interferer 
.its  i udependent ly  of  adjacent  interferers,  the  error  metric  will  be  equal 
in  magnitude  and  separated  by  A0.  The  most  remote  error  peak  (i.e.  the 

I irgest  | Y | location)  should  occur  for 

“i  ■ "t  1 (b“rbs  + <"»  - 1)  • f)  ■ 

where  BW  is  the  resolution  beamwidth.  Ideally,  these  standard  deviation 
K lib 

< -aroma  should  occur  for  A0  = 6 and  LL  = 40  at  the  positions  indicated  by 
toe  sin;  1 1 l arrows  in  figure  E.5.  It  is  seen  that  the  actual  extrema  loca- 

l ions  and  the  predicted  locations  (i.e.  the  arrows)  are  in  excellent  agree- 
ment lor  LL  = 40  and  A0  = 6°.  However,  A0  is  not  quite  large  enough  to 
produce  standard  deviation  of  error  peaks  of  equal  magnitudes. 

As  mentioned  earlier,  the  peak  SBT  bias  error  for  larger  array 
sizes  ( see  figures  E.5  through  E.10)  is  not  only  smaller  in  magnitude  but 
also  occurs  closer  to  the  true  target  bearing.  Unlike  the  standard 
d vial  ion  metrics,  it  is  quite  difficult  to  perceive  the  underlying  relation- 
siip  between  the  location  and  "shape"  of  the  peak  bias-error  and  the 
i .-soliit  ion  power  of  the  SBT  processor.  One  phenomena  worth  noting,  however, 
i i the  effect  of  very  widely  spaced  interferers  on  the  SBT  bias  error.  The 
i solul ion  beamwidth  of  the  40-hydrophone  array  is  small  enough  to  cause 
t lie  appearance  of  multiple  local  extrema  for  A0  = 6°  (see  figures  E.8  and 
I .10).  These  peaks  illustrate  the  "push-pull"  effect  that  directional 
noise  has  on  the  bearing  estimation  bias.  As  an  interferer  closes  on  the 
I irget  I rom  remote  bearings,  the  target  bearing  estimate  is  increasingly 
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"pul Icil"  or  biased,  in  the  direction  of  the  interferer.  When  the  point 
of  peak  bias  Is  reached,  the  interferer  is  effectively  "pushing"  the  bear- 
ing estimate  hack  toward  the  target.  In  figures  E.8  and  E JO  each  inter- 
fere?- in  the  clusters  of  2 and  4 interferers  are  seen  to  exert  a "push- 
pull"  i t le.ct  on  the  bearing  estimates  as  shown,  respectively,  by  the  two 
and  fom  local  maxima  (minima  for  Y < 0) . These  local  maxima  indicate  that 
at  some  Y-separatlon  each  interferer  acts  as  the  primary  cause  of 
csLfmal ion  bLas.  It  is  unclear  from  the  magnitudes  of  these  local  maxima 
to  Jus i what  extent  each  interferer  is  the  primary  cause  of  bias  error. 

Other  than  these  basic  qualitative  statements,  little  can  be  said  concerning 
the  relationship  between  peak  bias  error  location  and  the  array  resolution 
heamwi  d III. 

The  simulation  standard  deviation  metric  is  again  larger  than  the 
approximate,  theoretical  metric.  It  is  also  noted  that  there  is  some  dis- 
crepancy between  the  simulation  and  theoretical  bias  metrics  for  small  bias 
magnitudes.  In  the  previous  section,  a similar  discrepancy  occurred  at  low 
bias  magnitudes  but  appeared  to  decrease  in  severity  as  the  bias  magnitude 
Increased  (l.e.  figure  D.12,  SNR  = 2.0).  It  is  believed  that  these  dis- 
crepant ies  are  similar  — both  most  likely  resulting  from  the  simulation 
hearing  resolution  limit  and  the  small  (N  = 20)  sample  size. 

We  can  again  obtain  a better  idea  of  the  effects  of  extended  array 
size  on  hearing  estimation  error  by  examining  the  standard  deviation  of 
error  fur  remote,  peak  and  coincident  interferers  and  the  peak  bias  error 
.is  fmu  i ions  of  the  array  size.  The  standard  deviation  metrics  for  remotely 
positioned  interferers  is  shown  in  figure  6.9  as  a function  of  the  extended 
range  ol  array  sizes.  Since  the  error  metrics  appear  linear  in  logarithmic 


Kigure  6.9  - Standard  Deviation  of  Error;  Remotely 
Positioned  Interference 
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looid  I nates,  the  standard  deviation  metrics  are  exponential  functions  of 
1,1..  I'ur thermore,  the  slopes  of  the  theoretical  metrics  indicate  that  the 
standard  deviation  of  the  bearing  estimates  varies  approximately  in  inverse 
propoi  I ion  with  the  3/2  power  of  array  size,  LL.  On  the  other  hand,  the 
slope  of  the  simulation  metric  is  -1  indicating  a simple,  inverse  proportion- 
al ilv.  In  a previous  paper  [19],  MacDonald  and  Schultheiss  addressed  the 
Sill  hearing  estimation  performance  in  the  presence  of  an  incoherent  noise 

field.  They  indicate  that  in  isotropic  noise  the  standard  deviation  of 

-3/2 

bearing  error  J or  high  SNR  will  vary  as  a function  of  LL  [19,  p.  42, 
equation  38].  Consequently,  since  it  has  been  shown  that  remotely  posi- 
tioned iulcrferers  have  little  effect  on  the  estimation  standard  deviation 

(e.g.,  see  figure  5.5  and  the  relative  isotropic-noise  error  magnitudes 

-3/2 

in  f ij’iire  K.l),  we  must  conclude  that  LL  dependence  indicated  by  the 
theoretical  metrics  is  correct  and  that  the  simulation  metric  behavior  is 
in  i s I cud  i ng . 

Recall  that  there  was  an  abrupt  increase  in  the  simulation-target 
discrepancy  in  section  6.1  when  SNR  was  increased  from  1.0  to  2.0  (see 
figure  6.2).  The  behavior  of  the  simulation  metric  shown  in  figure  6.9 
reflects  a similar  phenomena,  i.e.  the  larger  arrays  provide  a higher  gain, 
or  rolieronl  signal  level,  and  as  before  cause  a larger  discrepancy.  The 
smaller  slope  of  Lhe  simulation  metric  is  attributed  to  a limit  on  bearing 
resolution  try  the  simulated  estimation  algorithm.  Note  that  there  is  a smaller 
sample  size  (20  runs)  and  fewer  beam  positions  per  estimate  (9)  for  LL  =*  20 
and  40  than  for  l,L  = 10. 
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i slope  of  -I  will  It-  at  LI.  = AO  this  slope  has  Increased  and  the  simulation 

iiu; trie  essentially  mirrors  che  theoretical  trend. 

The  theoretical  metrics  are  all  similar  in  their  LL-dependency , 

which  has  decreased  from  LL  ^ ^ that  we  observed  for  the  remote  inter- 

lerers  to  II.  ' ^ in  this  case  — a small  but  notable  change.  Note  also 

i hat  i lie  peak  error  is  dependent  on  number  and  separation  of  interferers. 

I'hls  contrasts  with  Jittle  or  no  dependence  on  MM  and  A0  shown  by  these 

metrics  for  remotely  positioned  interferers. 

For  the  case  of  coincident  target  and  interference,  the  standard 

deviation  of  error  metrics  are  shown  in  figure  6.11.  The  CR  bound  is 

again  seen  to  be  higher  than  the  theoretical  SBT  metric  for  some  LL  values 

and  a difference  in  the  slopes  of  the  two  SBT  metrics  is  again  apparent. 

A new  observation  is  the  tendency  of  the  standard  deviation  of  error 

metrics  to  increase  as  the  number  of  closely  spaced  interferers  increases. 

-3/2 

for  one  coincident  interferer,  the  error  is  proportional  to  LL  while 
lor  A "coincident"  interferers  It  is  approximately  proportional  to  LL  ^ . 

Thus,  there  is  a slower  improvement  in  estimation  accuracy  with  increasing 
array  size  for  larger  number  of  interferers. 

As  a final  indicator  of  the  effects  of  anisotropic  noise,  consider 
the  behavior  of  the  peak  SBT  bias  relative  to  array  size  as  summarized 
In  figure  6.12.  These  results  indicate  a large  improvement  in  bias  error 
r < > r larger  array  sizes.  A fourfold  increase  in  the  number  of  hydrophones 
shows  that  the  theoretical  bias  metrics  decrease  by  factors  of  between  3 
and  I while  the  simulation  metrics  decrease  by  factors  hetween  9 and  13. 

As  In  the  other  cases,  this  theoretical-simulation  bias  discrepancy  is 
at  tr llml  ed  to  the  sample  size  and  the  I.SE  line  fit. 
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Absolute  Peak  Bias  (degrees) 


CHAPTER  VII 


CONCLUSIONS 

/ . I - Summary 

In  many  source  bearing  estimation  problems,  highly  anisotropic 
noise  I ields  occur  as  a rule  rather  than  an  exception.  Such  problems 
generate  interest  in  the  performance  of  both  optimal  and  sub-optimal,  fixed- 
design,  sonar  array  processors.  Non-adaptive  optimal  processors  represent 
a design  goal  while  simplicity  and  economy  motivate  the  use  of  fixed-design, 
sub-optimal  processors.  Consequently,  this  research  has  investigated  the 
effects  of  highly  anisotropic  noise  fields  on  several  performance  metrics 
lor  the  maximum-likelihood  (non-adaptive)  estimator  and  the  conventional, 
split-beam  tracker.  This  study  has  addressed  target  tracking  problems 
where  the  signal-to-total-interf erence  ratio  is  large  (S/MI>1)  and  where 
true  target  bearings  are  approximately  broadside. 

We  have  considered  the  performance  of  the  optimal  bearing  estimator 
I rom  the  standpoint  of  the  maximum  likelihood  estimation  criterion.  The 
performance  of  the  maximum  likelihood  estimator  is  characterized  in  terms 
of  a bound  on  the  standard  deviation  of  estimation  error,  thus  avoiding 
problems  associated  with  solutions  to  the  maximum  likelihood  equation.  A 
particularly  pleasing  metric  for  target  tracking  problems  is  the  Cramer-Rao 
bound.  Chapter  2 sought  to  derive  a theoretical  expression  for  the  Cramer- 
Kao  bound  which  is  expressable  in  terms  of  the  inverse  noise  covariance 
matrix.  Although  the  algebraic  inversion  of  the  covariance  matrix  of 
highly  anisotropic  noise  is  intractable,  numerical  matrix  inversion  algo- 
rithms permit  evaluation  of  this  expression.  The  Cramer-Rao  bound  is  derived 
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for  a stationary,  Caussian  signal  and  noise  model.  All  propagating  compo- 
nents of  the  acoustic  field  are  assumed  to  be  plane  wave  fronts  arising 
from  sources  located  in  the  far-field  of  the  receiving  array.  The  noise 
field  consists  of  isotropic,  uncorrelated  ambient  noise  and  MM  directional 
inter I erences . 

The  effects  of  the  same  anisotropic  noise  field  on  bearing  estima- 
tion by  the  split-beam  tracker  (SBT)  is  examined  theoretically  in  chapter 
3.  The  sub-optimal  SBT  estimates  a source  bearing  from  the  cross-correlation 
lag  time  of  two  split  beam  outputs.  The  steering  angle  at  which  the  proc- 
essor output  crosses  zero  identifies  the  bearing  estimate.  Chapter  3 
develops  a theoretical  SBT  output  expression  and  derives  further  expressions 
for  the  mean  and  variance  of  the  SBT  output.  Under  restrictive  assumptions 
on  the  steering-angle  dependence  of  the  SBT  output  mean  and  variance,  an 
approximate  expression  for  the  bearing  estimation  variance  is  derived. 

In  the  absence  of  experimental  bearing  estimation  results,  a 
digitaL  SBT  simulation  is  developed  to  provide  statistical  measures  — the 
sample  mean  and  variance  — of  SBT  performance  . The  SBT  is  simulated  in 
the  frequency  domain  to  avoid  potential  problems  arising  from  digital 
models  of  certain  temporal  operations  (e.g.  time  differentiation).  This 
approach  is  facilitated  by  elimination  of  the  Fourier  transformer  require- 
ment through  direct  generation  of  the  Fourier  random  deviates.  Conse- 
quently, the  SBT  processor  algorithm  is  available  from  the  results  of 
chapter  3.  The  bearing  estimation  algorithm  is  configured  as  the  zero- 
crossing  of  the  LSE  line  fit  to  the  outputs  at  ISTEER  beam  positions 
around  the  target  bearing. 

Chapter  5 identifies  the  CR  bound  as  the  sole  performance  metric 
of  the  optimal  (Ml,)  estimator.  Also  defined  are  simulation  and  approximate, 
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theoretical  measures  of  SBT  estimation  bias  and  standard  deviation.  These 


optimal  and  sub-optimal  performance  metrics  are  evaluated  for  a baseline 
set  signal,  noise  and  array  parameters.  The  functional  dependence  of 
the  performance  measures  on  target-to- interference  separation  is  presented 
for  I,  2,  and  A interferers  with  mutual  interference  separations  of  3,  6, 
and  9 degrees.  For  a more  restricted  set  of  interference  parameter  values, 
the  performance  results  are  obtained  for  extended  SNR  and  array  size  para- 
meters. Chapter  6 presents  and  analyzes  the  performance  metrics  for  an 
SNR  of  0.5  and  2.0  (baseline  - 1.0)  and  subsequently  for  array  sizes  of 
20  anti  AO  hydrophones  (baseline  - 10  hydrophones) . 

The  major  contributions  of  this  research  are  described  below: 

1)  The  Cramer-Rao  bound  on  the  variance  of  any  unbiased  estimator  as 
derived  in  chapter  2 is  a new  contribution  to  the  literature.  This 
result  is  applicable  to  completely  arbitrary  (complex),  total-noise 
covariance  matrices  and  requires  only  the  assumptions  of  linear 
receiving  array  and  stationary,  Gaussian  processes.  It  is  thus  more 
general  and  widely  applicable  than  existing  results. 

2)  Although  it  is  a fairly  straightforward  procedure,  the  extension 
of  the  theoretical  and  simulation  SBT  error  expressions  in  chapters 

3 and  A to  Include  multiple  uni-directional  interferers  is  believed 
to  be  original  with  this  dissertation.  These  expressions  permit  further 
study  into  the  effects  of  changes  such  as  in  array  geometry  and  degrees 
of  anisotropy.  With  slight  modification,  it  should  also  be  easy  to 
address  arbitrary  (spatially  distributed)  interference  fields. 

3)  The  numerical  results  presented  in  chapters  5 and  6 are,  to  our 
knowledge,  the  first  quantitative  study  of  the  effects  of  multiple 


InLerl ecers  on  passive  sonar  bearing  estimation.  Previous  known  studies 
liave  been  restricted  to  singie  interferers  and  a special,  random- 
parameter  , mul ti-interf erer  case. 

7.2  - Conclusions 

The  numerical  results  presented  in  chapters  5 and  6 illustrate  that 
performance  levels  obtainable  in  the  presence  of  highly  anisotropic  noise 
fields  are  complex  functions  of  the  degree  of  noise  anisotropy.  The  degree 
of  noise  field  anisotropy  ranges  between  two  extremes.  At  one  extreme  is 
the  case  of  multiple  coincident  interferers  acting  as  a single  interference 
with  a level  of  MM- I.  The  other  extreme  is  MM,  widely  separated  inter- 
ferers each  affecting  the  estimation  performance  independently  (at  different 
hearings,  of  course).  The  transition  between  different  phenomena  at  each 
extreme  accounts  for  the  complex  behavior  of  the  bearing  estimation  per- 
formance metrics.  Specifically  we  can  make  the  following  conclusions: 
l)  All  three  standard  deviation  metrics  exhibit  similar  behavior  in 

each  of  three  distinct  regions  of  target-to-interf erence  separation. 

For  remotely  positioned  interference  clusters  (figure  5.5),  the  Cramer- 
Rao  bound  is  unaffected  by  the  directional  interference  thus  indicating 
a complete  Interference  rejection  by  the  optimal  processor.  The  SBT 
interprets  an  increase  in  the  number  of  remote  interferers  as  a slight 
increase  in  equivalent  isotropic  noise  level.  Again,  the  separation  of 
interferers  is  irrelevant.  At  intermediate  separations  of  the  target 
and  interference-cluster  a peak,  or  relative  maxima,  is  observed  for 
all  standard  deviation  metrics.  This  peak  occurs  when  the  target  and 
closest  interference  are  approximately  one  resolution  beamwidth  apart. 

The  degradation  of  the  SBT  performance  over  the  optimal  estimator  in 


this  Intermediate  region  is  primarily  affected  by  the  interference 
separation,  AO  (see  figure  5.6).  In  other  words,  the  performance  loss 
ol  l lie  SBT  relative  to  the  CR  bound  increases  dramatically  with  increases 
In  i lie  total  interference  power  received  from  bearings  approximately  one 
resolution  beamwidth  on  either  side  of  the  source  bearing.  When  the 
target  and  interference  cluster  are  coincident,  the  effect  of  decreasing 
interference  separation  is  again  dramatic.  For  small  separations,  the 
SBT  Interprets  an  interference  cluster  in  this  region  as  a component  of 
the  source  field.  The  effect  is  an  increase  in  the  received  coherent 
field  and  a corresponding  decrease  in  standard  deviation  of  error  (often 
below  the  CR  bound).  As  A0  increases,  the  amount  of  interference  power 
concentrated  one  beamwidth  to  either  side  of  the  source  Dearing  increases 
and  lienee  the  standard  deviation  of  error  increases  as  A9  increases 
(figure  5.7).  The  same  phenomena  affects  the  CR  bound  but  at  a slower 
rate. 

2)  Multiple  interfering  noise  sources  bias  the  source  bearing  estimates 

of  Hie  SBT  in  the  direction  of  the  interference  cluster.  Peak  bias  error 
occurs  when  the  closest  interference  is  within  one  resolution  beamwidth 
of  the  source.  Peak  SBT  bias  increases  linearly  with  decreasing  inter- 
ference separation  and  appears  to  be  directly  proportional  to  the  number 
of  interferers  (figure  5.11). 

J)  In  the  presence  of  multiple  interfering  noise  sources,  the  ML  and 

SBT  standard  deviation  metrics  are  approximately  inversely  proportional 
to  i lie  square-root  of  SNR  (figures  6.2  - 6.4).  For  isotropic  noise, 
the  same  (SNR)  dependence  occurs  [19]  and,  hence,  the  degree  of  noise 
anisotropy  is  seen  to  affect  only  the  overall  random  error  magnitude  and 
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not  its  functional  dependence  on  SNR.  The  peak  SBT  bias  error  is 
inversely  proportional  to  SNR  and  directly  proportional  to  the  number 
of  closely  spaced  interferers  (figure  6.5).  Consequently,  the  peak 
SBT  bias  is  approximately  inversely  proportional  to  the  signal-to-total- 
Interference  ratio,  S/MM* I,  arriving  within  a resolution  beamwidth  of 
the  target. 

'*)  The  ML  and  SBT  standard  deviations  of  error  for  multiple,  remotely 

posii toned  interferers  (like  the  error  for  isotropic  noise  [19])  vary 

Inversely  with  3/2  power  of  array  size  (figure  6.9).  The  peak  standard 

deviation  metrics  decrease  slower  with  increasing  array  size,  varying  in 
-1.37 

proportion  to  LL  (figure  6.10).  For  coincident  target  and  inter- 

-3/2 

ference  cluster,  the  array  size  dependence  ranges  from  LL  to  about 
- 1 

LL  allowing  a weakening  dependence  on  LL  as  the  number  of  interferers 
increases  (see  figure  6.11).  Also,  the  SBT  bias  metrics  exhibit  a 
dramatic  decrease  with  higher  resolution  (larger  LL)  — the  improvement 
being  most  significant  with  the  larger  numbers  of  interferers. 

>)  With  few  exceptions,  the  theoretical  and  simulation  measures  of 

SBT  performance  consistently  exhibit  similar  behavior  although  the 
simulation  standard  deviation  is  typically  higher  than  the  theoretical 
measure.  In  those  cases  where  a bias  metric  discrepancy  is  suggested, 
the  simulation  bias  tends  to  be  lower  than  its  theoretical  counterpart. 
Four  possibilities  are  suggested  in  section  5.3  as  contributing  causes 
of  this  discrepancy. 

7.3  - Recommendations  for  Further  Work 

A natural  extension  of  this  research  would  investigate  the  effects 
of  multiple,  high-level  interferers  and  distributed  (possibly  non-Gauss ian) , 


anisotropic  noise.  Since  it  is  apparent  that  frozen-design,  sub-optimal 
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processors  cannot  tolerate  nearby  high-level  interference,  nulling  or 
inter  I crence-cancel ling  sub-optimal  processors  and  adaptive,  maximum- 
likelihood  estimators  must  be  examined.  Although  considerable  research 
has  been  done  on  optimal,  adaptive  estimation,  we  are  aware  of  only 
MacDonald's  [18]  study  of  sub-optimal  processor  structures  in  the  presence 
of  a single,  high-level  interference.  Also  of  interest  is  the  possibility 
of  Improving  processor  discrimination  against  multiple,  coherent  inter- 
ferences (high  and  low  level)  with  three-dimensional,  planar,  and  non- 
uniform  linear  array  geometries.  Several  authors  (e.g.  [18],  [19],  [5]) 
have  examined  the  performance  of  non-uniform  linear  (usually  symmetric) 
pass i vi!  sonar  arrays  for  the  isotropic  and  single- interference  cases. 

The  methodology  of  this  investigation  could  also  be  expanded  to 
include  the  effects  of  multiple,  non-s tat ionary  interferers  and  non- 
stat  ionary  signals.  Propagating  components  are  often  non-stationary  due 
to  multipath  transmissions  and  platform  dynamics.  Noise  transients  and 
sonar  countermeasures  are  also  sources  of  non-stationary  interference. 

Such  an  investigation  might  prove  fruitful  for  certain  cases  of  non- 
staiionarity  (e.g.  non-stationary  noise  processes  describable  as  locally 
stationary  [3]). 

The  Cramer-Rao  bound  was  derived  in  chapter  2 under  the  tacit 
assumption  that  the  source  bearing  was  the  only  unknown  parameter  in  the 
sign.)  I -plus-noise  model.  This  is,  of  course,  untrue  in  practice.  It 
would  be  interesting  to  examine  the  behavior  of  the  CR  bound  for  the  case 
of  multiple,  unknown  parameters.  This  is  especially  so  in  light  of  the 
sub-optimal  processor's  assumed  ignorance  of  the  directional  noise  para- 
meters. Such  an  investigation  would  require  derivations  (such  as  in 
chapter  2)  of  the  J2  elements  of  the  Fisher's  information  matrix  [33] 
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where  .1  is  the  number  of  unknown  signal,  noise  and  interference  parameters. 
The  <l<  bound  on  the  variance  of  any  unbiased  estimator  of  the  jth  unknown 
par a mol or  is,  by  definition,  the  ( J j ) element  of  the  inverse  Fisher's  in- 
formation matrix. 

Another  possible  area  of  investigation  might  address  active  sonar 
processor  performance  in  the  presence  of  multiple  interferences.  Here, 
source  range  and  range-rate  are  also  unknown,  but  required  parameters. 
Although  ti>e  directional  noise  sources  are  commonly  signal  dependent  in 
tills  case  (target-like  scatterers),  sonar  countermeasures  possibly  con- 
stitute signal-independent  interference.  Background  noise  transients  may 
also  represent  directional,  signal- independent  interference. 


APPENDIX  A 


DERIVATION  OF  THE  LOG-LIKELIHOOD  RATIO 

Equation  (2-9)  defines  the  Fourier  characterization  of  the  total, 
received  data  vector  as 

Xjj  - (vJu.q),  X2(jo,0).  ....  XL(Jo*0),  XjtJujj).  ....  Xjju^)]  ' 

(A-l) 

A T “Jw.t 

where  X^jWj)  - / x^(t)  e dt  (A-2) 

0 

and  x^( t)  is  the  real  waveform  received  by  the  £th  hydrophone  over  the 

-A 

Interval  [0,  T].  The  definition  in  (A-2)  of  the  components  of  Xp  as 
linear  transformations  of  the  real,  Gaussian  temporal  processes  suggests 
that  X^(ju>i)  is  also  a Gaussian  random  process.  The  components  of  X^  are 
by  definition,  however,  complex  random  processes.  The  formal  definition 
of  a complex,  Gaussian  process  [6]  requires  that  (1)  the  real  and 
imaginary  parts  of  the  complex  process  be  jointly  Gaussian  and  (2)  the 
covariance  matrix  of  the  real  and  imaginary  parts  be  expressable  in  an 
anti-symmetric,  block-matrix  form.  By  the  definition  (A-2),  the  real  and 
imaginary  parts  of  X (jw^)  are  jointly,  Gaussian  when  x^(t)  Is  a Gaussian 
process.  The  second  requirement  also  follows  from  the  definitions  in 
(A-l)  and  (A-2).  The  interested  reader  is  referred  to  chapter  2 of  the 
dissertation  of  U.  Bangs  [2]  for  a detailed  proof  of  the  validity  of  this 
requirement.  Specifically,  Bangs  proves  that 
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e |r«  [4]  • in  [4']|-  e j to  ['4]  • r«  [4]} 

where  Re  [•]  and  Im  [■]  denote  the  real  and  imaginary  parts  and  where 
is  defined  in  (A-l). 

Since  the  data  vector  satisfies  the  requirements  of  a complex, 
vector  Gaussian  random  process,  its  pdf  can  be  shown  to  be  [6] 


«(4|\)=  ("|\| ) 


exp  [ -Xjj  Xy 

( P 


where  is  assumed  to  have  a zero-mean, 

p = 0 or  1, 

Hj  denotes  the  signal  plus  noise  hypothesis, 

denotes  the  noise  only  hypothesis, 

is  the  covariance  matrix  of  X^  under  H and, 
P P 

is  the  determinant  of  . 

P P 

The  covariance  matrices  are  defined  as 


(A-3) 


v ■ eI44*|»p| 


(A-4) 


where,  from  (A-l)  and  (A-2),  an  arbitrary  element  of  is 

P 


E |X.«VXn(K>!  • / / dt  l d"  E |xm 


-j(m  t - iii.il) 

x (t)x  (u)  e 

m n ' (A-5) 


Since  x(t)  Ls  assumed  to  be  zero  mean  and  stationary  we  can  define  the  cross- 
correlation between  phones  as 


R (x)  = R (t-u)  * E | x (t)x  (u)  j 
mn  inn  | m n ) 


and  rewrite  (A-5)  as 


( ) T -ju>,T  T 

E ! X.<J“1>  j ■ ‘ 


~JWiT  T -j(“i  - uk)u 


du 

(A-6) 


When  T is  much  greater  than  the  correlation  time  of  x (t)  then  the  second 

m 

integral  in  (A-6)  can  be  approximated  as  T times  the  Kronecker  delta 
1 1 i = k 

S±k  ' |o  1 O 

and  the  first  integral  approximately  equals  the  cross-power  spectral 
density  between  the  m and  nth  phones.  Hence,  we  car  write 


E X (jui  ) xWja.  ) = T X (u>.)6,, 

J m i n k l mn  i ik 


(A-6) 


where  X^^w^)  is  the  cross-power  spectral  density  at 

Equation  (A-6)  allows  us  to  express  the  L x N covariance  matrix 

K„  in  the  following  block-diagonal  form 
P 


4 <"o> 

p 

= 0 
P 


<V. 

P 


(A-7) 


P 


where  0 is  an  L x L null  matrix. 


Ky  (u>.)  = E j X(uk)  X (a)i)  j is  the  L x L cross-covariance  matrix 


at  the  ith  frequency,  and 


X(w1)  is  the  vector  of  L components  of  X^  corresponding  to 
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The  block  diagonal  form  of  the  total  covariance  matrix  in  (A-7)  is  the 
objective  of  the  Fourier  characterization  of  the  received  data  given  by 
equations  (A—  1)  and  (A-2). 

With  Xp  rewritten  in  the  form 


Xp  “ j^X  , X (to^),  X (w^)  > •••»  X j 


(A-8) 


and  Kp  in  the  block  diagonal  form  of  (A-7),  it  is  easy  to  show  that  the 
P 

probabilities  defined  in  (A-3)  become 


V / N 

f 

\_l  ( 

N ‘ (’i 

P 

“Pl 

' v-1 


i-1 


(A-9) 


The  log- likelihood  ratio  is  by  definition, 

pfclv  9t) 


- tn 


•(\|Ho) 

and  upon  substitution  of  (A-9)  becomes 


kkw  J:  iSjVH 
' I ' v-1  \\\) 


(“i^”1  Mmp  - X *(»1)^  X (mi)  J 


N 

■ l ln^ 

i=l 


0|  * _lt^  -k  * -1 

U “X  Kj,  1 X + X K„  X 

Kp  I "^1  ^0 


(A-10) 


where,  for  simplicity,  the  frequency  subscripts  have  been  suppressed  and  the 
"prime"  notation  used  in  defining  the  submatrices  in  (A-7)  has  been  dropped. 

Evaluation  of  the  inverse  of  is  facilitated  by  the  following 

matrix  identity  for  matrices  of  the  form 


L 
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A 


b b 


+ C : 


-1  -1 
£ b b * £ 


* 

1 + b 


C b 


(A-n) 


Equation  (A-ll)  is  easily  verified  after  a minor  amount  oi  algebraic 
manipulation.  Recall  from  equation  (2-20)  that  the  data  covariance  under 
hypothesis  is 

JSn  = TS  a a +1^  (A-12a) 

and,  when  no  signal  is  present, 


(A-12b) 


where  is  the  covariance  matrix  (at  the  implied  frequency)  of  the  Fourier 
characterization  of  the  total  received  noise  vector.  Combining  (A-ll)  and 


(A-12a), 


1 


= 


-1 


TS 


a a 


-1 

h 


1 + TS  a 


h 


-1  a_».  -1 

4 a * Sn 


(A-13) 


where 


2 


o 


T 


A TS 

= i r*— -i 

1 + TS  a K a 


Substituting  (A-12b)  and  (A-13)  into  (A-10) , yields 


- j, 


£n 


IM 


2 * _1  -v  _i  * -1  _i 

°T  X % 3 a X 


(A- 1*0 
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The  ratio  of  matrix  determinants  In  (A-14)  is  evaluated  by  using 
the  matrix  identity 


and  equation  (A-13): 


I \l 

r^7[  ■ is  i ■ i*»  i 

I | i-l  2 -i  -V  _ i * _1 
|^n|  ' " °T  a a 

It  2 -*  -i  * -1  , 

’ |i  ' °T  a “ 4 | 

-1 

The  m,  n - element  of  the  matrix  a a K..  is 


-J.  -i  * 
a a 


ra,n 


L -j  (m  - p)  A 
1 e 
p=l 


m,n 


= e~j(»A)  £ e+j(p4) 

p=l 


1,  2, 


fA~15) 


L 


(A-16) 


, d 

where  A = u>k  - cos  0 , k = 0,  1 N and 

K^150  is  the  p,  n-element  of  1 . 

Since  the  row- index  m in  (A-16)  appears  only  as  an  exponent,  the.  deter- 
minant in  (A-15)  can  be  further  reduced,  since  the  addition  of  a constant 
times  the  rth  row  of  any  matrix  to  the  sth  row  of  that  matrix  does  not 

change  the  value  of  its  determinant.  Hence,  if  each  element  of  the  first 

2 -1  -l  * -1  -1(r-l)A 

row  of  _I  - oT  a a is  multiplied  by  -e  J and  added  to  the 

corresponding  elements  in  the  rth  row  for  r = 2,  3,  . . . , L,  the  following 

expression  is  obtained: 
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Expansion  of  the  determinant  down  the  1st  column  yields 


V 


V 


1 - l e-J(s*>C 

s=l 

2 L L 

1 - °T  l I e 

s=l  p=l 


-jsA  +jpA  ps 
% 


- 1 - ^ l I '+3pl  ></s 

s=l  p=l 


2 _»  * _l  _i 

= 1 - oT  a a 


--  r 

TS  |_ 


, * -1-1 
1 + TSa  iSfj  a-TSa 


‘ * -1  -»  1 
h 3J 


T 

TS 


(A-19) 


Substitution  of  (A-19)  in*-o  (A-14)  yields  the  desired  expression  for  the 
log-likelihood  ratio 


ftl'r)  = j,  t"|(S)*4*V‘ 


-V  _i*  -1  _v) 

a a Kjj  Xj 


(A-20) 


where 


°T  = 


TS 


-A  * - l-l 

1 + TS  a a 


S = S(u)j)  is  the  signal  power  spectral  density  at 
-1 

is  the  inverse  total-noise  covariance  matrix  at  id^ 

T is  the  observation  interval, 

-l 

a is  the  target  delay  vector  at  and 

X is  the  L- dimensional  received-data  vector  also  at  the  frequency 
Equation  (A-20)  appears  as  equation  (10)  in  chapter  2 of  the  text. 
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APPENDIX  B 


DERIVATION  OF  EQUATION  (2-21) 


The  CR  bound  requires  the  expected  value  of  the  second  derivative 
of  the  log-likelihood  ratio  at  the  arbitrary  frequency  i within  the  band- 
width  of  the  observation  vector  random  process,  X (u>,): 


*n !x_  +!_ 

3y  TS  Bp 


V * -1  -i  _i  * -1  _i 

A a ^ X X a 


— 1 ft  ft  ”1  \ _>  ^ “]_j  dO_  ^ ^ — ] | t 

+ 2—  • A K„  XX  ^ a + 2-E,  ^ X X K A 

3p  N 

2 * * * — 1— x _i  it  — 1—1  2 —1  _j 

+ °T  a AAI^XX  ^ X,,  XXt^AA. 


„ 2 -J  * * -1  ->  -i  * -1  _» 

+ 2 oT  a A XX  ^ A a 


(B-l) 


From  equal- ion  (20)  of  chapter  2 the  expectation  of  the  matrix  X • X is 


I — 1 * | -i  -i  * 

E | X X j = TS  a a + 


The  expectation  of  the  second  term  in  (B-l)  is 


(B-2) 


5T  _»  * -1  \-S.  _i  *1  -1  _»  d °T  ->  * -1  r — i -Vt  -1 

^ 3 K*  E X X a . a ^ I TS  a a + ^ a J 


*|  -1  _»  3 °T  -j  * -1 

K»,  a *•  — a K,, 


2 2 2 2 

3 °T  -»  * “!-»-»*  ~1  —i  3 °T  _i  * ~1  _s 

= TS  — 2“  a a a a + — z1  a a 

3m  3p 


2 2 

3 °T  _»*  -1_> 
3p 


6) 


(B-3) 


The  last  term  of  (B-l),  upon  taking  its  expectation,  becomes 
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2 * * -1 
2 oT  a 


* -1  f *)  -1 

A ^ E | X X \h  A 


2 L * * -i  -*  -*  * -1  -x  * -l 

2 °T  a A (TS  a a ^ A a + a A A a 


2 TS  o„ 


-x  * * -1  _x 


a A 


A** 


2 -x  * * -1  -V 


+ 2oTa  AK^  Aa 


(B-6) 


Substituting  (B-3),  (B-4) , (B-5),  and  (B-6)  Into  (B-l)  yields,  for  the 


expected  value  of  the  second  derivative  of  t 


i-lil  _ * 


2.  2 2 
3 


I 3p  ) 3p 


o 


In  M-  + 


3p 


+ 2 


+ (TS) 


/3°t\  TS  * f * -1  _1  I _» 

[77/7/  l4  *■  + 1H* 

[-x  * * * -!_»_>*  -1 

3 A A Kg  a + a AA  aj 


2 X * * -1  -»  2 2 -.  * * -1 


+ 2 TS  o,|%  A ^ 


+ 2ot*  A ^ A a (b-7) 


The  underlined  terms  of  equation  (B-7)  can  be  combined  for  further  slmpll- 

f icatlon,; 


[_»***  “1  _x  _x  * “1  _x  I 2 _x  * * “1  _» 

■ A*?.  a + a is  AAaj+2cT.  A a 

2 ( _»  * * -i  _»  / _»,  * -l  _x\[_i  * * * -i  _»  _» * -i  — n 

“°x)2a  A An  Aa  + (1  + TSa  a)la  A A An  a+a  A A a 

/-'*  -i-»\  r-x  * -i  -x-i 

’ (a  An  a 1 Re  [a  An  AA  a j 

2 ( -i  * * r * -i  -i  i -x  jv*r  * -l  -i-i  _x) 

°Tja  A [A  An  +An  Aj  a + a [A  An  aJ  AaJ 


2ot  (TS) 


CB-8) 


Finally,  substituting  (B-8)  into  (B-7)  yields  the  desired  result 
equation  (2-21)  in  chapter  2. 
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APPENDIX  C 

DERIVATION  OF  CLOSED  FORM  EXPRESSION 
FOR  DOUBLE  SUMMATION  IN  EQUATION  (4-28) 


The  first  double-summation  in  equation  (4-28)  is  quickly  evaluated 


as, 


L 

2 L 


I Ii  «-p)  ■ 5 I p 

1=1  pJ±+l  2 £-i  2 L , 

2 P=»2+l 


L L 

2 2 


I t-l 


1=1  r=l 


L 

2 


_ L ? L 
2 *■  ' 2 


L 

8 


(C-l) 


r=l 


A similar,  but  longer,  process  is  followed  for  the  second  double  summation 
in  (4-28): 

L L 


l 


l l . a - p>  cos  [a  - p)nl  - * i l,  a-  P>Lja  - p)!V  - p)!1j 

*1  L J l=i  p=^+-i  ( ) 

k k L L 1 

2 2 / -ja-r-f)fl 


~ik  n( 


S e J2  "|J  U+W  J e^10  l re’3'0-  \ | J e'^J 


+ % e+J^lj  Uita  l e+3t!1  -I  l re+3t0  - i J J .+3ri!  j 

U r * r 2 £ r 


j sin 


I ^e_:}£n  J e+jrn  - j.  e"jto  J re+Jrnj  - \ 


cos  2 n 


UO  v -jr« 


e- 


(C-2) 


^ L 


j 
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where  all  summation  indices  run  from  1 to  y • 
Due  to  the  equality 


t. 

I iJte  = _«  _JL_  7 
l-i  J an  £ 


CC-3) 


there  is  only  one  unique  summation  in  (C-2) , for  which  the  following  Identity 
applies: 

Un  _ j(Hl)  ,i„  (|  • f ) 


1 - e 


sin  | 


(C-4) 


where  we  define. 


4k  • i: 

sin  (§) 


(C-5) 


Substitution  of  (C-4)  and  (C-3)  into  the  bracketed  term  in  (C-2)  yields. 


h (<T)ta  l eJO>  - £ e'J&  Ii 


-H1)  «] 


“1 

r il 

+ 2)n/ 

-i(L  + 2)n 

H-  c 

J 

: 

> 2 rl 

ac  2 r 

Jan  e 

/L  4- 

i, 

/ +1 

[L  + 2) 

-j(— 

( . ac  J 
Man  e 

' 4 ’iH1) 

2 j c 
J an 


From  (C-5),  it  is  easy  to  obtain 


(C-6) 
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APPENDIX  D 


NUMERICAL  RESULTS  FOR 
EXTENDED  SNR  VALUES 

This  appendix  presents  the  bearing  estimation  error  results 
obtained  for  extensions  of  the  baseline  SNR  to  0.5  and  2.0.  The  three 
standard  deviation  of  error  metrics  appear  in  figures  D.l  - D.7  and  the 
two  SBT  bias  error  metrics  appear  in  figures  D.8  - D.14  as  functions  of 
the  target- to- interference  separation  Y,  where 

Y • (H)  «(£)  (cos  «t  - co’  ei ) * 

= mean  interference  bearing, 

0^  = target  bearing  90°, 

T = 0.128  seconds,  KK  =*  13, 
d = 20  feet,  c = 5000  ft/sec, 

LL  = 10  hydrophones,  INR  * 0.1. 

In  figures  D.l  - D.7  the  following  legend  applies: 

CR  Bound  on  the  Standard  Deviation  of  Error  of 
Any  Unbiased  Bearing  Estimator, 

Approximate  Theoretical  SBT  Standard  Deviation  of 
Error, 

Simulated  SBT  Standard  Deviation  of  Error  and 
90%  Confidence  Belts  (20  replications  per  sample) 
for  SNR  - 0.5  and  SNR  - 2.0,  respectively. 


I 

The  solid  line  shown  in  the  bias  error  results  of  figures  D.8  - 
n.14  represents  the  approximate  theoretical  SBT  bias  metric  while  the 

! 

simulated  SBT  bias  results  are  depicted  as  above  for  the  random  errors. 

i 

i 


standar  V 

DE.VI4T/OH 


SNR  » 0.5 


SNR  * 1.0 


* 


/t.Af 


SNR  - OS 


S}4R  * t o 


— j— 

i.o 


Figure  D.5  - Standard  Deviation  of  Error  Results  for  Extended 

SNR  Values;  MM  - 4,  A0  - 3* 
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’ ‘ ‘ I 

Figure  D. 6 - Standard  Deviation  of  Error  Results  for  Extended  I 

SNR  Values;  MM  - 4,  A0  - 6* 
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Figure  D.14  - SBT  Bias  Error  Results  for  Extended  SNR  Values 
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APPENDIX  E 

NUMERICAL  RESULTS  FOR 
EXTENDED  ARRAY  SIZES 

This  appendix  presents  the  bearing  estimation  error  results 
obtained  for  extensions  of  the  baseline  array  size  (LL)  to  20  and  40 
hydrophones.  The  three  standard  deviation  of  error  metrics  appear  in 
figures  E.l  - E.5  and  the  two  SBT  bias  error  metrics  appear  in  figures 
E.6  - E.10.  The  appropriate  baseline  values  of  the  metrics  are  repeated 
in  figures  E.l  - E.10  for  comparison.  The  random  and  bias  errors  are 
presented  as  functions  of  the  target-to-interference  separation,  Y,  where 

' ' *■  T-ei)’ 

■ mean  interference  bearing, 

0T  =*  target  bearing  = 90° 

T “ 0.128  seconds,  KK  ■ 13, 
d *=  20  feet,  c * 5000  ft/second 
SNR  * 1.0,  INR  »0.1  . 

In  figures  E.l  - E.5  the  following  legend  applies: 

CR  Bound  on  the  Standard  Deviation  of  Error  of 

Any  Unbiased  Bearing  Estimator, 

_ Approximate  Theoretical  SBT  Standard  Deviation  of 
Error, 

Simulated  SBT  Standard  Deviation  of  Error  and 
•fc  A and  0 90%  Confidence  Belts  for  LL  ■ 10,  20,  and  40 

hydrophones,  respectively.  For  figures  E.l  - E.10 
the  SBT  simulation  sample  sizes  are: 


JJ  J 


11  1 
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CONFIGURATION 


SAMPLE  SIZE 


LL  - 

10, 

MM  - 

1 

60 

LL  - 

10, 

MM  » 

2,  4 

40 

LL  - 

20, 

40 

20 

The  solid  line  shown  In  the  bias  error  results  of  figures  E.6  - E.10 
represents  the  approximate  theoretical  SBT  bias  metric  while  the  simu- 
lated SBT  bias  results  are  depicted  as  above  for  the  random  errors. 


Figure  E.2 


Standard  Deviation  of  Error  Results  for  Extended 
Array  Sizes;  MM  - 2,  AO  - 3* 
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